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ABSTRACT 


A  software  defined  radio  is  a  much  more  flexible  platform  than  traditional, 
hardware  implemented  radios.  By  implementing  radio  functions  in  software,  and  putting 
those  functions  on  a  Field  Programmable  Gate  Array  (FPGA)  chip,  users  will  have  the 
ability  to  download  mission  specific  radio  capabilities.  This  thesis  examines  a 
fundamental  piece  of  the  receiver,  the  Phase-Lock  Loop  (PLL),  simulates  a  software 
PLL,  and  investigates  the  effects  of  fixed-point  versus  floating  point  mathematics 
required  for  an  FPGA  based  PLL.  With  a  fixed-point  PLL  simulator,  figures  of  merit 
such  as  lock-time,  lock  range,  and  pull-in  range  are  determined  for  typical  signal-to-noise 
ratio  (SNR)  levels. 
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EXECUTIVE  SUMMARY 


The  explosion  of  wireless  and  PCS  serviees  over  the  last  deeade  has  created 
numerous  incompatible  air  interface  standards.  A  subscriber  to  one  service  will  find  his 
phone  rendered  useless  when  roaming  in  the  coverage  area  of  another  service,  using  a 
different  standard.  These  competing  transmission  formats  each  might  have  their  own 
unique  modulation  type,  multiple  access  technique,  error  control  methods,  call  set-up  and 
handoff  protocol  and  voice  compression  algorithms.  The  military  has  seen  similar 
interoperability  of  radio  standards  between  tactical  radios  used  by  coalition  forces  in 
Desert  Storm. 

The  need  for  a  flexible  communication  platform,  capable  of  interfacing  with  the 
numerous  standards  and  formats  has  become  apparent.  A  software  radio,  implementing 
traditional  radio  functions  in  software,  gives  the  user  this  flexibility.  The  ideal  software 
radio  shown  in  Figure  (1)  would  digitize  the  entire  received  signal  spectrum  using  a  high 
speed  Analog-to-Digital  Converter  (ADC),  perform  all  demodulation,  data  protocol  and 
processing  functions  using  a  general-purpose  digital  signal  processor  (DSP). 


Figure  1.  Ideal  Software  Defined  Radio.  (From  Ref  [3].) 
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Due  to  limitations  on  A/D  conversion  speed,  sampling  the  direct  RF  spectrum  is 
not  typically  an  option.  A  practical  software  radio  will  usually  incorporate  an  RF  front 
end,  which  filters  and  downconverts  a  portion  of  the  received  signal  spectrum  to  much 
lower  IF  frequency  such  as  21.4,  70  or  160  MHz.  These  frequencies  can  be  digitized 
directly  with  current  state  of  the  art  A/D  converters.  This  technique  is  called  “IF 
Sampling”. 

The  digital  signal  processing  (DSP)  functions  of  the  software  radio  can  be 
implemented  either  with  general-purpose  DSP  chips  or  with  reconfigurable  Field 
Programmable  Gate  Arrays  (FPGAs).  In  either  case,  the  processing  functions  are  in  the 
form  of  software,  available  for  download  to  the  DSP  engine.  The  resulting  unit  would 
have  the  capability  to  be  reconfigured  for  any  radio  signal  format.  This  gives  the  user  the 
ability  to  download  mission  specific  radio  requirements,  using  the  same  platform  for 
numerous  radio  applications. 

This  thesis  will  focus  on  the  implementation  issues  for  an  important  signal 
processing  function  common  to  most  communications  receivers;  that  is,  the  Phase-Lock 
Loop  (PLL).  The  PLL  can  take  several  forms  such  as  the  Costas  Loop  for  carrier 
recovery  and  tracking,  the  early-late  gate  synchronizer  for  baud  timing  recovery,  and  the 
delay-locked  loop  (DLL)  for  spreading  sequence  tracking  in  spread  spectrum  systems.  In 
each  case,  the  same  PLL  loop  theory  presented  in  this  thesis  applies. 

In  short,  the  PLL  is  a  feedback  loop  device,  which  locks  onto  a  received  signal, 
meaning  it  synchronizes  its  output  in-phase  and  frequency  with  its  input.  The  PLL  can  be 
broken  down  into  its  three  component  parts:  1.)  the  phase  detector  (PD),  2.)  a  loop 
lowpass  filter,  and  3.)  a  voltage  controlled  oscillator  (VCO)  or  numerically  controlled 
oscillator  (NCO),  the  latter  being  used  in  a  software  PLL  version.  Figure  2. shows  a 
simple  PLL  schematically. 
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Figure  2.  Simple  Phase-Loek  Loop 

While  a  traditional  analog  reeeiver  implements  the  PLL  components  exclusively 
in  hardware,  a  software-defined  radio  requires  these  components  to  be  available  as 
downloadable  software,  using  a  digitized  signal  as  the  input.  While  a  software  simulation 
of  a  PLL  has  been  a  reality  for  some  time,  the  advent  of  fast  Field  Programmable  Gate 
Array  (FPGA)  technology  makes  this  a  useful  concept  for  high  data  rate  signals  (those 
greater  than  roughly  2  Mbit/sec  data  rate).  Currently,  the  software  components  can  be 
programmed  onto  an  FPGA  and  used  real-time  for  DSP  functions. 

When  deciding  on  DSP  implementations,  one  must  consider  whether  fixed-point 
or  floating  point  arithmetic  and  number  representation  will  be  used.  The  benefit  of 
floating  point  is  the  large  dynamic  range  associated  with  the  floating-point  number 
representations.  When  implementing  algorithms  in  floating  point,  the  designer  typically 
does  not  have  to  worry  about  issues  such  as  rounding  or  truncation  error,  or  numeric 
overflow.  The  disadvantage  of  floating  point  is  the  increased  computational  resources 
required  and  the  processing  speed  limitations.  For  this  reason,  floating  point 
implementations  of  receiver  algorithms  on  general  purpose  DSP  are  limited  to  relatively 
low  data  rate  signals. 

One  of  the  reasons  FPGAs  are  fast  enough  to  be  used  is  the  fact  that  they  use 
fixed-point  representations  of  numbers  rather  than  floating  point.  Fixed-point 
representation  is  a  much  more  efficient  way  for  a  computer  to  do  arithmetic,  because  its 
essentially  a  binary  representation  of  a  decimal  number.  Hence,  a  computer  can  do 
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arithmetic  at  a  much  faster  rate  using  fixed-point  numbers.  With  the  computationally 
heavy  nature  of  any  DSP  application,  the  efficiencies  of  fixed-point  arithmetic  add  up  to 
significant  savings  in  time.  The  trade  off  is  that  representing  numbers  in  this  way  either 
reduces  the  range  of  the  number,  or  its  precision.  As  an  example,  using  8-bit  fixed-point 
numbers,  a  range  of -256  to  256  only  has  a  precision  of  2.0.  For  a  precision  of  .01,  the 
range  of  numbers  available  are  -2  to  2.  This  constraint  needs  to  be  realized  when 
programming  using  fixed-point  arithmetic. 

The  aim  of  this  research  is  to  model  a  software  phase-lock  loop,  observe  its 
performance,  and  convert  it  into  a  fixed-point  implementation  to  determine  the  effects  on 
important  performance  such  as  lock  time  and  pull-in  range.  The  questions  that  required 
answering  were 

•  How  many  bits  are  needed  in  a  fixed-point  implementation  for  acceptable 
performance? 

•  What  type  of  errors  do  a  fixed-point  implementation  introduce  into  the 
output  signal? 

•  What  effect  does  fixed-point  arithmetic  have  on  figures  of  merit  of  a  PLL? 

•  What  kind  of  signal-to-noise  ratios  (SNR)  are  required  to  lock  the  PLL  in 
an  acceptable  time 

The  source  code  used  initially  was  MATLAB,  but  to  convert  the  model  into  fixed- 
point,  the  MATLAB  extension  Simulink  was  used.  The  final  simulation  model  used  16- 
bit  fixed-point  arithmetic  and  locked  with  acceptable  SNR’s,  with  only  small  errors  to  the 
output  signal.  Lock  time  and  ranges  were  not  changed,  however  a  smaller  pull-in  range 
than  the  floating-point  equivalent  was  encountered.  Errors  at  the  output  of  the  PLL  due 
to  the  fixed-point  implementation  are  observed,  but  nothing  of  critical  size.  The  fixed- 
point  model  performed  comparable  to  an  analog  or  floating-point  model  in  all  ways,  and 
could  used  as  a  basis  to  build  a  PLL  on  a  FPGA  or  DSP  chip. 
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I.  INTRODUCTION 


A,  BACKGROUND 

The  revolution  in  wireless  technology  in  the  last  decade  has  created  a  need  for  a  fast, 
flexible,  and  light  radio,  compatible  with  the  numerous  transmission  formats  and  standards.  The 
concept  of  a  software  radio  has  been  developed  to  fill  that  need.  The  ideal  software  radio 
digitizes  the  radio  spectrum  at  the  receive  antenna,  providing  all  demodulation,  decryption,  and 
signal  processing  in  software.  Such  a  tool  would  have  the  capability  to  download  appropriate 
software,  depending  on  the  mission  and  the  transmission  format  required.  It  could  be 
reconfigured  to  accommodate  any  RF-band  modulation  or  data  format  or  transmission  standard, 
capable  of  operating  within  any  communications  network. 

To  realize  the  ideal  software  radio,  the  entire  signal  spectrum  would  be  digitized  at  the 
antenna  by  an  analog-to-digital  converter  (ADC).  A  practical  SW  radio  will  usually  incorporate 
an  RF  front  end  which  filters  and  downconverts  a  portion  of  the  received  signal  spectrum  to  much 
lower  IF  frequency  which  can  be  digitized  directly  with  current  state  of  the  art  A/D  converters. 
The  digital  signal  processing  (DSP)  is  then  performed  in  software  in  a  reprogrammable  Field 
Programmable  Gate  Array  (FPGA).  Finally,  the  demodulated  signal  is  sent  through  a  digital-to- 
analog  converter  (DAC)  to  generate  the  audio  or  video  output  if  required.  The  ideal  software 
radio  architecture  is  shown  in  Figure  1,  where  the  DSP  functions  are  done  on  an  FPGA  chip.  The 
radio  concept  described  above  could  be  reprogrammed  to  accommodate  any  radio  standard  or  air 
interface  by  downloading  the  appropriate  software  algorithm. 

The  signal  processing  functions  required  by  the  software  radio  includes  tuning,  fdtering, 
demodulation  and  decryption.  Tuning  is  accomplished  by  mixing  the  digitized  signal  with  a 
digital  local  oscillator  to  down-convert  the  signal  to  baseband.  This  mixer  is  simply  a  signed 
multiplier,  sample  by  sample  in  the  digital  implementation.  Typically,  a  complex  baseband  signal 
representation  is  used.  The  mixer  will  output  the  in-phase  and  quadrature-phase  components  in 
sine  and  cosine  waveforms.  This  mixing  process  creates  unwanted  frequencies,  specifically  at 
twice  the  receive  frequency.  This  spectral  component  is  removed  using  a  finite  impulse  response 
filter  (FIR).  The  digital  filtering  also  bandlimits  the  samples  to  the  bandwidth  of  the  signal  of 
interest  (a  process  called  “channelization”). 
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Demodulation  is  performed  by  a  variety  of  methods,  depending  on  the  modulation  type. 
This  research  assumes  demodulation  will  be  performed  by  phase  locking  to  the  incoming  signal. 
An  error  signal  is  fed  back  to  the  tuner  to  adjust  the  local  oscillator  frequency  in  a  closed  loop 
system.  This  feedback  system  is  simply  a  digital  or  software  phase-lock  loop  (PLL).  The 
implementation  of  such  a  PLL  is  the  focus  of  this  research. 

1,  The  Phase-Lock  Loop 

A  simple  phase-lock  loop  is  pictured  in  Figure  2  and  can  be  seen  to  consist  of  three 
component  parts:  1.)  the  phase  detector,  2.)  a  loop  fdter,  and  3.)  a  voltage  controlled  oscillator 
(VCO)  or  numerically  controlled  oscillator  (NCO).  Its  purpose  is  to  lock  on  to  the  frequency  and 
phase  of  the  input  signal.  In  this  case,  the  purpose  of  the  PLL  is  to  create  a  phase  coherent  as  a 
local  oscillator  in  the  receiver  for  demodulation. 

The  phase  detector  or  phase  comparator  compares  the  phase  between  the  input  signal  and 
the  output  signal.  It  generates  a  signal  proportional  to  the  phase  error,  or  difference  between  the 
two  signals  phase.  This  can  be  done  numerous  ways;  the  simplest  is  to  multiply  the  two  signals 
together. 

The  loop  filter  block  is  a  low-pass  filter  that  removes  the  high  frequency  terms  the  come 
from  the  multiplication  of  the  input  and  output  signal,  leaving  only  the  phase  error. 

The  NCO  takes  the  phase  error  from  the  loop  filter  output  and  adjusts  its  output  sinusoidal 
signal  to  force  the  error  to  zero.  This  adjusted  signal  is  the  feedback  signal  that  goes  into  the 
phase  detector,  producing  a  second  phase  error  and  the  process  repeats. 

This  type  of  PLL  is  a  second-order  system.  The  error  eventually  settles  to  zero,  but  the 
output  is  a  damped  oscillation.  This  oscillation  is  governed  by  the  parameters  of  the  PLL. 


2.  Fixed-Point  Arithmetic 

Measurements  of  physical  quantities  can  take  on  many  numerical  representations.  For 
example,  the  number  one  thousand  can  be  represented  by  1000,  1E3,  10  ,  or  one  thousand.  In  this 

case,  the  same  quantity  is  represented  using  four  different  syntaxes.  Another  example  of  the  same 
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quantity,  represented  in  different  scales,  is  the  boiling  point  of  water.  Water  boils  at  100  degrees 
Celsius,  212  degrees  Fahrenheit,  373  degrees  Kelvin,  or  671.4  degrees  Rankine. 

Determining  an  appropriate  seale  and  representation  depends  on  many  faetors.  Suppose 
you  need  to  measure  the  speed  of  a  vehiele.  The  numerieal  values  have  a  limited  range  for  this 
applieation.  The  slowest  it  ean  travel  is  0  mph  at  a  dead  stop  and  the  top  speed  of  the  vehiele  has 
been  determined  to  be  150  mph.  If  an  8-bit  unsigned  integer  is  required,  values  in  the  range  of  0- 
255  are  possible. 

A  typieal  approaeh  would  be  to  assign  one  bit  per  mile-per-hour,  making  the  integer  0  a 
dead  stop  and  150  the  vehiele’s  top  speed.  This  seheme,  while  easy  to  eonvert,  negleets  the  use  of 
the  integers  151-255,  wasting  40%  of  the  number  range. 

Another  approaeh  would  be  to  set  the  integer  zero  to  a  dead  stop  and  the  integer  255  to  the 
top  speed.  This  seale  gives  us  mueh  greater  preeision,  0.58823  mph  per  bit,  beeause  all  256 
values  are  used  in  the  vehiele  speed.  The  eonversion,  however,  requires  a  division  of  1.7,  a 
relatively  expensive  operation  for  fixed-point  proeessors. 

The  trend  of  reeent  teehnology  is  to  implement  eontrol  systems  and  digital  signal 
proeessing  funetions  on  digital  hardware.  Within  digital  hardware,  numbers  are  represented  as 
either  floating-point  or  fixed-point  data  types.  The  number  of  bits  used  to  represent  both  data 
types  is  a  fixed  word  size.  The  range  of  fixed-point  representation  is  mueh  smaller  than  for 
floating  point,  thus  to  avoid  overflow  and  quantization  errors,  fixed-point  representation  must  be 
sealed.  If  floating-point  numbers  ean  effectively  represent  real  world  values  with  virtually  no 
error,  why  use  a  fixed-point  based  implementation?  The  answer,  of  course,  is  cost,  size  and 
processing  speed. 

A  fixed-point  hardware  platform  is  architecturally  much  simpler  than  its  floating-point 
counterpart.  This  means  cheaper  manufacture  of  the  product.  In  addition  to  manufacturing 
savings,  if  scaled  properly,  fixed-point  arithmetic  can  be  significantly  faster,  saving  computation 
time. 

Because  the  logic  circuits  of  fixed-point  hardware  are  less  complicated  than  that  of 
floating  point,  the  chip  size  can  be  much  smaller,  reducing  power  consumption.  This  means 
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smaller  batteries,  and  reduced  heat  as  a  by-product,  removing  the  need  for  an  expensive  and  bulky 
heat  sink. 

For  a  software  radio,  programmed  on  a  FPGA  where  speed  is  a  critical  limitation,  and  size 
and  battery  life  are  a  concern,  a  fixed-point  implementation  is  the  only  plausible  solution. 


B,  OBJECTIVES 

This  research  focused  on  the  effects  of  fixed-point  arithmetic  on  a  simulated  PLL, 
specifically,  what  number  of  bits  are  required  to  successfully  implement  a  PLL,  how  does  this 
affect  the  figures  of  merit  of  a  PLL  such  as  pull-in  range,  lock  time,  and  lock  range,  and  to 
determine  the  signal-to-noise  ratios  (SNR)  needed  for  a  fixed-point  PLL  to  be  effective  in  a 
software  radio. 

These  objectives  were  accomplished  using  MATLAB’s  Simulink  software.  A  PLL 
simulation  using  only  fixed-point  arithmetic  was  implemented,  taking  as  input  a  noisy  signal. 
This  simulation  was  done  using  16-bit  fixed-point  arithmetic.  Simulation  results  were  consistent 
with  theoretical  results,  with  the  fixed-point  implementation  simply  adding  “noise”  to  the  signal, 
causing  the  output  signal  to  have  small  errors  due  to  the  quantization  effects  of  fixed-point 
arithmetic. 

C.  RELATED  RESEARCH 

The  related  research  is  in  identifying  the  feasibility  and  requirements  for  other 
components  to  make  a  software  radio  realizable.  Some  of  this  research  is  being  done  concurrently 
with  this  research.  The  first  is  a  digital  delay  lock  loop  (DDLL);  the  second  is  investigating 
digital  downconversion  and  channelization. 

1.  Digital  Delay  Lock  Loop 

For  a  fully  functional  software  radio,  a  real-time,  spread-spectrum,  signal-processing 
block  is  desired.  An  integral  function  of  a  direct  sequence  spread  spectrum  receiver  is  recovering 
the  underlying  narrowband  data  through  a  procedure  called  dispreading.  The  key  component  of 
dispreading  is  a  DDLL.  A  DLL  generates  an  exact  replica  of  the  spreading  sequence  generated  at 
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the  transmitter,  allowing  the  reeovery  of  the  narrowband  data.  This  researeh  was  being  eondueted 
by  Captain  Samuel  Laboy,  USMC. 

2,  Digital  Downconversion  and  Channelization 

A  basie  building  bloek  of  any  all-digital  reeeiver  is  the  downcoversion  of  a  real  signal  at 
an  intermediate  frequeney  to  eomplex  in-phase  and  quadrature-phase  eomponents  at  baseband. 
The  downcoversion  component  of  the  software  radio  consists  of  a  complex  NCO,  local  oscillators 
that  act  as  mixers,  and  digital  filter.  Investigation  in  the  feasibility  of  digital  downconversion  and 
performance  characteristic  compared  to  an  analog  equivalent  is  being  pursued  by  Lieutenant 
Michael  Snelling,  USN.  The  aim  of  this  research  is  to  build  a  working  MATLAB  model  and  to 
analyze  the  results  of  the  simulation. 


D,  THESIS  ORGANIZATION 

This  thesis  is  organized  to  mirror  the  research.  Chapter  2  is  an  overview  of  the  operation 
of  the  PLL  and  derives  the  figures  of  merit  for  analysis  of  the  final  model.  Chapter  3  is  a  tutorial 
on  fixed-point  arithmetic  and  explains  the  advantages  that  it  offers  over  floating-point  arithmetic. 
Chapter  4  is  the  culmination  of  this  research.  It  begins  by  describing  the  building  of  a  floating¬ 
point  PLL  and  then  the  procedures  for  converting  it  to  a  fixed-point  model  in  Simulink.  Analysis 
of  the  performance  of  the  fixed-point  PLL  according  to  the  derived  figures  of  merit  derived  in 
Chapter  2  is  also  done  in  Chapter  4.  The  conclusions  for  the  feasibility  of  a  FPGA  based  PLL  and 
recommendations  for  future  research  are  in  Chapter  5.  The  Appendix  gives  the  MATLAB  code 
and  Simulink  model  for  the  floating-point  and  fixed-point  phase  lock  loops. 


5 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


6 


II.  PHASE  LOCK  LOOP  BASICS 


A,  PRINCIPLES  OF  THE  PHASE-LOCK  LOOP 

A  phase-lock  loop  is  a  circuit,  or  software,  designed  to  track  a  given  Reference  signal  in 
both  frequency  and  phase.  Its  applications  are  far  reaching  from  AM/FM  radio  demodulation,  to 
television  sets,  to  coders  and  decoders.  The  PLL  has  the  ability  to  synchronize  its  output  in-phase 
and  frequency  with  an  input  signal,  meaning  the  phase  error  between  the  PLL’s  output  and  the 
input  signal  is  zero,  or  remains  constant.  If  a  phase  error  is  introduced,  by  a  phase  change  or  a 
frequency  change,  the  PLL’s  feedback  control  mechanism  adjusts  the  oscillator’s  output  to 
account  for  it. 

To  get  a  basic  understanding  of  the  PLL,  the  principles  of  the  linear  PLL  will  be 
examined,  which  is  shown  schematically  in  Figure  3.  It  consists  of  three  blocks:  the  phase 
detector  (PD),  the  loop  filter  (LF),  and  the  VCO/NCO.  For  ease  of  Reference,  the  same  notation 
and  signal  names  used  by  [I]  will  be  used  here.  Consequently  the  frequency  of  Uy{t)  in 
radians/second  is  ft),,  (O^  is  the  frequency  of  WjCOand  the  phase  error  or  difference  in- 

phase  of  Wj(t)  and  WjCO  .  With  signals  defined,  a  closer  look  at  each  functional  block  is  in  order. 


Figure  3.  Phase-Lock  Loop  with  Signals  of  Interest 
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1.  Phase  Detector  (PD) 

The  PD  takes  as  its  input  two  signals,  the  reference  signal,  Wj(t),  and  the  output  of  the 
VCO,  ^2(0  •  This  block  simply  compares  the  phases  of  these  two  signals  and  produces  a  signal 
(t)  proportional  to  the  phase  error  0^ ,  specifically, 

uAt)  =  KA  (2.1) 

where  is  the  PD  gain  in  volts.  Obtaining  the  signal  Uj{t)  can  be  as  simple  as  a  multiplier. 

The  PD  adopted  in  the  simulation  uses  the  in-phase  and  quadrature-phase  portions  of  the  input 
and  output  signal.  To  see  how  this  works,  lets  assume  a  phase  difference  between  the  u^{t)  and 
(t)  of  0^ .  Defining  the  in-phase  and  quadrature-phase  components  as  /;  =  [/;  cos(<Wit) , 

=  f/j  sin(<Wjt) ,  /j  =  t/2  cos(<W2^  +  Qi  -  ^2  sin(<W2^  +  ,  so  that 

Mj  {t)  =  f/j  (cos(<Wjt)  -I-  i  sin(<Wjt)) 

M2  (0  =  U2  (cos(<W2^  +  )  +  *  sin(<W2t  +  )) 

where  t/j  and  t/j  1^^  amplitude  of  the  respective  signals  and  is  related  to  the  PD  gain  . 
Using  the  trigonometric  identities 

cos  X cos  y  =  ^[cos(.r  -I-  y)  -I-  cos(.r  -  y)]  andsin  vsin  y  =  ^[cos(v  -  y)  -  cos(.r -I-  y)] 

and  a  little  algebra,  0^  can  be  extracted.  For  ease  of  computation,  we  assume  0)^  =  0)2.  If  this  is 

not  the  case,  the  PD  output  has  an  additional  frequency  error,  but  this  ac  component  will  be 
filtered  out  by  the  loop  filter  explained  below. 

I  =  I,^l2+Q,^Q2 

=  t/j  *t/2[cos(tyjt)cos(ft>2^  +  6’J  +  sin(c9jt)sin(<»2^  +  ^J] 

1  (2.2) 

=  —U^*U2[cos(((O^+a)2)t  +  0J  +  cos(0J  +  cos(0J-cos(((O^+(O2)t  +  0J] 

=  U,*U2COS(0J. 

Similarly, 
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Using  equations  (2.2)  and  (2.3), 


=  U,*[/,sin(0J. 


(2.3) 


0^  =  tan  ' 


I  1) 


As  desired,  the  PD  output  signal  (t)  is  the  phase  error. 


(2.4) 


2,  Loop  Filter 

As  alluded  to  earlier,  if  the  frequeneies  of  the  two  signals  differ  or  alternate  PD 
implementations  are  used,  the  output  signal,  of  the  PD  will  have  an  unwanted  ac 

component.  This  ac  signal  is  superimposed  on  the  desired  dc  component  representing  the  phase 
error.  To  remove  the  ac  component,  a  simple  low  pass  loop  fdter  is  used. 

To  implement  a  loop  filter,  numerous  strategies  exist.  While  a  high-order  finite  impulse 
response  (FIR)  filter  will  ensure  the  ac  component  is  removed,  the  delay  in  response  is  too  great, 
where  locking  onto  a  signal  as  fast  as  possible  is  needed.  A  first-order  low-pass  filter  has  a  quick 
response  and  successfully  removes  the  unwanted  oscillations.  The  most  common  loop  filter  used 
is  called  an  active  PI  filter  (PI  =  proportional  +  integral,  taken  from  control  theory  and  named  due 
to  the  fact  that  it  has  a  pole  at  s  =  0,  hence  acts  as  an  integrator)  [1].  Taken  from  the 
corresponding  RC  circuit  filters  and  implemented  in  software  using  their  transfer  functions,  the 
filter  has  the  Bode  plot  depicted  in  Figure  4  [1].  The  transfer  function  for  the  PI  filter  is  given  by 

F(.)  =  l±^a  (2.5) 

5Tj 

where  in  the  analog  circuit  world,  xi  and  X2  are  RC  time  constants  of  the  circuit  fdter  and  are 
determined  by  the  loop  bandwidth  of  the  system,  which  is  in  turn  a  function  of  the  noise. 


9 


|F| 


T^2 


Figure  4.  Bode  Plot  of  Active  PI  Filter.  (From  Ref  [1]). 

3.  Voltage  Controlled  Oscillator  /  Numerically  Controlled  Oscillator 

The  NCO  takes  as  its  input  the  output  from  the  loop  filter,  Uj-  (t) ,  which  is  proportional  to 
the  phase  error.  By  using  an  appropriate  scaling,  or  gain  constant  ,  of  the  NCO,  we  can  adjust 
the  output,  MjCO?  to  account  for  the  phase  error  in  the  two  signals.  The  NCO  has  a  center 
frequency.  This  is  the  oscillator  frequency  without  any  adjustments.  This  value  is  selected  to 
reflect  approximate  working  frequencies  of  the  system.  The  PLL  will  properly  lock  onto  a  signal 
within  a  percentage  range  of  the  center  frequency.  Call  this  frequency  roo-  An  adjustment  made 
on  the  frequency  of  the  NCO  will  be  an  offset  from  wq-  Thus  the  output  of  the  NCO  has  a 
frequency, 

ft>2  (0  =  <x>o  +  (0-  (2-6) 

The  discrete  nature  of  the  software  system  requires  the  output  of  the  NCO  to  be  a  phase, 
as  opposed  to  frequency,  and  thus  by  definition  the  phase  is  given  by  the  integral  over  the 
frequency  variation, 

=  (2.7) 

Using  a  look-up  table,  the  output  of  the  system  is  taken  as 

=cos(6>2(0) 

^li^QuadraturePhase  ~  Sin(02  (0) 


10 


Either  the  in-phase  or  the  quadrature-phase  ean  be  used  as  the  oseillator,  but  with  the 
implementation  deseribed  above,  both  are  needed  as  inputs  to  the  PD. 


B.  OPERATION  OF  THE  PHASE  LOCK  LOOP 

To  understand  the  operation  of  the  PLL  and  how  the  bloeks  work  together,  assume  first 
that  the  system  is  initially  loeked,  i.e.,  0)^=  0)2,  and  examine  three  types  of  signals  at  the  input;  a 
phase  step,  a  frequeney  step,  and  a  frequeney  ramp.  Taken  from  [1],  the  signals  are  shown  in 
Figure  5. 


M\.f 

V  vy 

.1  -a^  D  D5  1 


Figures.  Exeiting  Funetions  as  Applied  to  Input  of  a  PFF.  (a)  Phase  Step, 
(b)  Frequeney  Step,  (e)  Frequeney  Ramp.  (From  Ref  [1]). 


The  aim  is  to  eharaeterize  these  signals  in  terms  of  their  phase,  and  proeeed  with  analysis 
using  a  generie  phase  term  in  the  signal.  We  assume  an  input  signal  of  a  sinusoid  as  previously 
done,  so  that  =  +  ,  where  is  a  generie  phase  signal  that  represents  the 

various  exeitation  input  signals. 

For  the  phase  step,  whieh  represents  a  phase-modulated  signal,  simply  performs  a 

step  ehange  at  t  =  0,  and  is  given  by, 
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6>i(0  =  A^>a(0 


(2.9) 


where  0  <  A^>  <  2k  is  the  magnitude  of  the  phase  ehange  and  u{t)  is  the  step  funetion. 

The  next  ease  is  an  example  of  a  frequeney-modulated  signal.  This  signal  has  a  step 
ehange  in  frequeney  at  t  =  0 ,  whieh  is  given  by  the  inerement  A<w .  Thus  the  eorresponding  input 
signal  is  given  by 

M,  (t)  =  sin((<w,  +  AQ)tu{t))) 

=  sin(<W]t  +  Acotu{t)). 

Thus 

e,{t)  =  Am.  (2.10) 

This  phase  signal  is  just  a  ramp  funetion  for  t  >  0 . 

The  final  signal  is  one  whose  frequeney  inereases  linearly  with  time.  As  shown  in  Figure 
5,  the  frequeney  is  a  ramp  funetion,  and  the  rate  of  inerease  in  frequeney  is  the  slope  of  the  line. 
If  the  slope  of  the  ramp  funetion  for  frequeney  is  Am,  then  the  total  frequeney  of  Mj(t)is 
<w,  +  Amt .  By  definition,  the  frequeney  of  a  signal  is  the  first  derivative  of  its  phase  with  respeet 
to  time, 

m,  +Am  =  — 
dt 

Thus  the  phase  of  the  signal  at  time  t  is  the  integral  of  its  angular  frequeney  over  the  time  interval 
0  <  T  <  t ,  and  so  the  input  signal  (t)  can  be  written  as 


M,  (t)  =  t/j  sin  j(<»j  +  Amt)dt 

0 


=  t/j  sin 


( 

md-t-Am 

V 


Consequently,  the  phase  signal  is  given  by 


ddt)  =  Am— 


(2.11) 
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1.  The  Transfer  Functioni 


It  is  often  helpful  in  examining  the  operation  of  any  system  to  derive  and  understand  the 
transfer  funetion,  H{s),  of  the  system.  This  will  be  the  approaeh  in  understanding  the  dynamie 
operation  of  the  PLL.  The  transfer  funetion  is  defined  as  the  Laplaee  transform  of  the  output 
divided  by  the  Laplaee  transform  of  the  input.  For  a  PLL,  we  want  to  relate  the  phase  signals  of 
the  output  and  input,  thus 

=  ^  (2.12) 

0,(s) 

where  0,  (5)  is  the  Laplaee  transform  of  the  phase  signal  6.{t)  .  To  build  this  mathematieal  model 
of  the  system,  assume  initially  the  PLL  is  loeked,  so  that  <w,  =  <^2  •  This  being  the  ease,  the  output 
of  the  PD  is  the  phase  error  6^  =  6^  times  the  PD  gain  eonstant  ,  i.e.. 


whieh  implies 


^d(t)  =  KA 


(2.13) 


(2.14)2 


The  next  bloek  in  the  PLL  is  the  loop  filter.  The  transfer  funetion  of  the  PI  aetive  filter 
was  already  given  in  Equation  (2.5).  The  last  bloek  to  derive  the  transfer  funetion  is  the 
VCO/NCO.  Reeall  that  the  NCO  adjusted  the  output  frequeney  of  u^it) , depending  on  its  input 
from  the  loop  fdter.  A  negative  input  redueed  the  output  frequeney  and  a  positive  input  inereased 
it.  The  eenter  frequeney  of  the  NCO,  0)^ ,  was  the  starting  point  for  this  adjustment.  The  output 

frequeney,  (O^it) ,  was  defined  by 


ft>2(0  =  ^0+'^oW/(0 


(2.15) 


where  K^Uf  (t)  is  the  variation  in  the  frequeney.  However  we  want  the  output  phase  of  the  signal, 

not  its  angular  frequeney.  Using  the  definition  of  the  phase  as  the  integral  over  the  frequeney 
variation, 

1  Transfer  function  analysis  and  derivation  done  using  [1]. 

2  Note  here  we  assume  over  a  small  time  period,  the  phase  error  does  not  change  dramatically  and  is  thus 
approximated  to  be  a  constant. 
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(2.16) 


^2(0  =  Kq  ^Uf{t)dt, 

The  Laplace  transform  of  an  integral  is  1/s,  thus 

=  (,)  ^  H,,,is)  =  ^.  (2.17) 

s  s 

The  closed-loop  block  diagram  of  the  transfer  function  is  shown  in  Figure  6. 


Figure  6.  Block  Diagram  of  PLL  Transfer  Function 


From  control  theory,  the  closed  loop  transfer  function  for  the  entire  system  is  given  by 


_  K,K,F{s) 

\  +  H,^{s)F{s)H,,^{s)  s  +  K,K,F{sy 


(2.18) 


Simplifying  this  to  standard  notation.  Equation  (2.18)  becomes 


H{s)  = 


^  °  (1-1- 572) 

F 


2  , 

5’  -\-s 


'kaF 


-I- 


(2.19) 


V 


J 


^2 


Control  theorists  put  this  in  terms  of  natural  frequency,  and  damping  factor,  ^ ,  in 

order  to  use  a  standard  equation  which  can  easily  be  understood  and  studied.  The  equivalent 
transfer  function  equation  in  normalized  form  is  given  by 


H{s)  = 


2sCo)„+0)l 

s^+2sCo)„+0)y 


(2.20) 
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where 


^  2 

To  ealeulate  the  time  response  of  the  transfer  funetion,  the  inverse  Laplaee  transform  of 
(2.20)  is  required.  The  time  domain  equivalent  equation  beeomes 


whieh  is  simply  a  damped  oseillation. 

To  get  a  feel  for  the  dynamie  operation  of  the  PLL  for  the  three  Referenee  inputs 
deseribed  above,  we  look  at  the  transient  response  of  the  transfer  funetion  in  normalized  form  for 
different  values  of  (O^  and  ^ .  Then  finding  a  desired  response,  values  to  build  the  PLL  ean  be 
determined. 

a.  Phase  Step  Applied  to  Input 

To  see  a  typieal  transient  response,  a  phase  step  was  applied  to  the  input  at  t  =  0 , 
eorresponding  to  the  first  type  of  input  above,  so  that 

6>i(0  =  A^>m(0  . 

Figure  7  shows  the  error  time  response  of  the  transfer  funetion  versus  time  for  various  values  of 
^  .3  A  natural  frequeney  of  200  rad/see  is  used.  The  value  of  <W„only  seales  the  x-axis,  not  the 

3  The  time  response  can  be  obtained  by  taking  the  inverse  Laplace  transform  of  the  transfer  function.  For  this 
transfer  function  the  time  response  is  given  by  Equation  (2.21) 
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response  of  the  system  or  its  oharaeteristies.  To  note  is  the  fact  that  as  the  error 

approached  zero.  This  can  also  be  derived  from  the  transfer  function  by  finding  the  error  transfer 
function,  defined  by 


liXs)  = 


0,(5) 


S 

s  +  K,K,Fis) 


(2.22) 


and  using  the  final  value  theorem  of  the  Laplace  transform  which  states 


0^(00)  =  lim^^o  50^(5)  =  O. 


Figure  7.  Time  Response  of  PLL  for  different  values  of  ^  and  =  200  . 


b.  Frequency  Step  Applied  to  Input. 

Again  applying  the  final  value  theorem,  a  frequency  step  corresponds  to  a  transfer 
function  input  of 

0,(5)  =  Aco!  s 
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Performing  the  same  ealeulation  on  the  phase  error  transfer  funetion, 


0^(00)  =  lim^^o 


Ao) 

+  2sCco„+co^„ 


Aco 


(o: 


(2.23) 


This  value  approaehes  zero  for  small  frequeney  steps  and  high  gain  loops,  whieh 
inereases  the  value  of  (O^ 

c.  Frequency  Ramp  Applied  to  Input 

Using  similar  teehniques  for  a  frequeney  ramp,  the  results  are  less  eneouraging. 

The  error  is  9^°°)  =  Ad)l  co\ ,  where  Am  is  the  rate  of  ehange  of  the  input  signal  frequeney.  Thus 

for  a  large  rate  of  ehange  in  frequeney,  the  PLL  unloeks.  Experimen-tation  with  PLL’s  has 
shown  that  a  praetieal  design  limit  to  the  rate  of  ehange  in  frequeney  is, 

=-^.  (2.24) 


C.  PERFORMANCE  MEASURES  OF  THE  PLL 

For  a  full  derivation  of  the  following  key  PLL  parameters,  one  is  eneouraged  to  see 
Referenee  [1].  For  the  purposes  of  this  researeh,  it  is  suffieient  to  define  the  parameters  and  state 
the  results  for  a  given  PLL,  and  shed  some  light  on  the  eompeting  parameters  and  their 
interaetion.  The  foeus  will  be  on  live  key  parameters  that  govern  the  dynamie  performanee  of  the 
PLL 

•  The  loek  range, 

•  The  loek  time,  Tl 

•  The  pull-in  range,  A(Op 

•  The  pull-in  time,  Tp 

•  The  pull-out  range  AtOp^  . 

Additionally,  a  figure  of  merit  for  the  PLL  of  loop  bandwidth.  Bp  will  be  given 
eonsideration  during  the  study  of  PLL  performanee  in  the  presenee  of  noise. 
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1.  The  Lock  Range 

The  lock  range  is  defined  as  the  range  of  frequency  offset  between  the  center  frequency  of 
the  VCO  and  the  reference  frequency  in  which  the  PLL  locks  within  one  single-beat  note  between 
reference  and  output  frequencies.  We  assume  that  the  PLL  is  initially  unlocked  and  switched  on 
at  time  t  =  0 .  To  derive  this  figure  of  merit,  we  assume  a  simple  PD  of  multiplying  the  input 
reference  signal  and  the  NCO  output  signal  to  produce  the  phase  error.  If  this  is  the  case,  then  for 
a  frequency  offset  of  A<y , 

{t)  =  sin(A<yt)  -l-  higher-frequency  terms  (2.25) 

where  the  higher-frequency  terms  are  discarded  due  to  the  loop  filter.  Looking  at  the  output  of 
the  loop  fdter,  the  result, 

Uf  (t)  =  I  T’(A<y)  I  sm(A<»)  (2.26) 

is  obtained.  The  output  of  the  LF,  Uj  (t) ,  is  just  an  ac  signal  causing  a  frequency  modulation  of 

the  NCO,  with  a  peak  frequency  deviation  of  \  F{Aq))  \ .  For  lock  to  occur,  this  peak 

frequency  deviation  must  be  less  than  the  frequency  offset  Ao) .  If  the  offset  is  larger  than  this 
frequency  deviation,  lock  cannot  occur  in  a  single  cycle.  The  lock  range  can  be  determined  by 
determining  when  this  peak  frequency  deviation  is  just  as  large  as  the  frequency  offset.  Thus  the 
equation  for  locking  becomes 

Aw,=K,K,\F(AwJ\.  (2.27) 

This  is  a  non-linear  equation,  but  to  solve  it  a  practical  approximation  can  be  made.  From  the 
Bode  plot  of  the  loop  filter  transfer  function  in  Figure  4,  the  lock  range  is  greater  than  the  cut-off 
frequencies  I/Tj  and  l/fj.  Thus  a  conservative  approximation  is  |  F(A<y^)  |=  Tj  /  .  Using 
Equation  (2.20),  the  lock  range  becomes 

Ao)^  -  2C(jO„ .  (2.28) 

Simulations  and  experiments  have  shown  that  this  is  a  conservative  approximation  and  can  be 
used  confidently  in  the  design  process.  Clearly  a  larger  lock  range  is  desirable. 
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1.  Lock  Time 


Lock  time  is  the  time  the  PLL  takes  to  lock  onto  a  signal  when  it  is  initially  unlocked, 
given  that  the  reference  signal  is  within  the  lock  range.  When  ^  <  1 ,  the  transients  of  the  damped 
oscillation  die  out  after  one  cycle,  thus  the  lock  time  ean  be  approximated  accurately  as 


(2.29) 


2,  The  Pull-in  Range 

The  lock  range  is  a  subset  of  the  pull-in  range,  whieh  is  defined  as  the  range  of  frequeney 
offset  in  which  the  PLL  will  eventually  lock  after  a  number  of  cycles.  This  is  dependant  on  the 
type  of  loop  filter  used.  For  the  aetive  PI  filter,  the  pull-in  range  is  infinite,  i.e., 

>oo,  (2.30) 

thus  any  referenee  frequeney  input  to  the  PLL  will  eventually  be  locked  onto.  For  different  types 
of  loop  fdters,  with  a  finite  gain  at  de,  the  pull-in  range  is  deereased.  See  [1]  for  further  details  on 
pull-in  range  for  different  filter  types. 


3.  Pull-in  Time 

Pull-in  time  is  the  time  required  for  a  PLL  to  look  onto  a  frequeney  within  the  pull-in 
range.  Defining  Aa)^  as  the  frequeney  offset  a)^-a)2,  the  pull-in  time  is  determined  to  be 


AcDq 

16  C(^l 


(2.31) 


Note  that,  as  Ao)g  o° ,  the  pull-in  time  approaches  infinity.  The  larger  the  offset,  the 
longer  the  pull-in  time. 


4,  The  Pull-out  Range 

The  pull-out  range  is  defined  as  that  frequeney  step  which  causes  the  PLL  to  unlock  if 
applied  to  the  reference  input.  Direot  oaloulations  of  the  pull-out  range  are  not  possible,  but  in 
[5],  oomputer  simulations  have  arrived  at  the  following  approximation 
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AO)po='^.SQ)SC  +  d- 


(2.32) 


The  relationship  of  the  frequeney  ranges  is  depicted  in  Figure  8.  Starting  from  the  center 
frequency  of  the  system,  the  lock-range  is  the  smallest  of  the  range  parameters.  Next  is  the  pull¬ 
out  range,  and  finally  the  pull-in  range  has  the  largest.  The  range  of  values  outside  the  pull-in 
range  is  called  dynamically  unstable  due  to  the  fact  when  operating  dynamically,  the  system  will 
not  lock  if  a  reference  frequency  in  that  range  is  encountered.  Inside  the  pull-in  range,  and 
outside  the  lock  range  is  called  conditionally  stable.  This  means  the  system  will  eventually  lock 
on  frequencies  in  this  range,  but  could  unlock  temporarily  with  a  reference  frequency  in  this 
range. 


Figure  8.  Relationship  of  frequency  ranges  of  a  PLL.  (From  Ref.  [1].) 
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D,  PERFORMANCE  OF  THE  PEL  IN  NOISE 


An  in-depth  analysis  of  the  performanee  of  the  PLL  in  noise  is  outside  the  seope  of  this 
researeh;  however,  some  insight  into  what  ean  be  expeeted  from  a  PLL  operating  in  the  real  world 
is  required.  The  following  definitions  are  used  and  assumptions  made  in  this  diseussion  of  noise. 

1.  SNRi=\0\og{Pg/P^),  where  P^  is  the  referenee  signal  power  and  P^  is  the  noise 
power,  is  the  signal-to-noise  ratio  at  the  input  of  the  PLL. 

2.  All  noise  is  assumed  to  be  Additive  White  Gaussian  Noise  (AWGN),  meaning  that  it 
has  a  fiat  power  speetral  density. 

3.  A  low-pass  pre-filter  of  bandwidth  B.  is  implemented  prior  to  the  PLL. 


4.  The  noise  bandwidth  is  defined  as  =  ^\H{jco)f  Jtu,  where //(jtu)  is  the  transfer 

0 

funetion  of  the  PLL.  Inserting  equation  (2.20)  for  H{jco) ,  the  integral  solves  to 


Bl=  — 
2 


f 


r  +  — 


(2.33). 


5.  The  signal-to-noise  ratio  at  the  output  is  defined  by 

=  Sm,  i .  (2.34)4 

IB^ 

This  figure  is  helpful  in  determining  how  often  the  PLL  unloeks  due  to  noise,  and  the 
SNR^  required  to  ensure  a  loek. 


Reeall  that  the  loek  range,  A<w^,  from  Equation  (2.28)  is  proportional  to  the  natural 
frequeney  and  is  desired  to  be  as  large  as  possible;  henee,  theoretieally,  designing  the  PLL  with  a 
<w„as  large  as  required  is  an  easy  task,  but  with  the  introduetion  of  noise,  and  Equation  (2.33),  the 
larger  ,  the  larger  is  5^ .  This  suggests  a  trade-off  between  redueing  the  noise  suffieiently  and 
a  larger  loek  range,  whieh  is  exaetly  the  ease.  To  inerease  the  loek  range  of  the  system,  a  larger 
SNR  must  be  used.  If  only  a  small  SNR  is  available,  5^  must  be  made  small,  whieh  reduees 
A(0^. 


4  For  a  derivation  of  SNR;,  one  is  direeted  to  Ref  [1]. 
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The  same  is  true  of  the  loek  time.  This  parameter  is  inversely  proportional  to  co^ ,  and  thus 
larger  values  of  (O^  reduee  .  Again,  by  redueing  is  inereased. 

Praetical  experiments  with  PLLs  have  shown  that  for  stable  operation,  we  need  an 
SNRi^  >  6  dB  [1],  A  system  eould  eventually  loek  with  a  lower  SNRj^ ,  but  the  phase  jitter  would 
eause  frequent  unloeks  and  the  system  would  be  of  little  value.  This  leads  to  the  final  important 
figure,  that  of  how  often  a  PLL  system,  on  average,  will  unloek.  is  defined  as  the  average 

time  interval  between  two  unloeks  of  the  system.  Figure  9  is  taken  from  [2],  which  depicts  as 
a  function  of  SNR^.  For  high  SNR^,  very  difficult  to  find  or  measure,  thus  Figure  9  is 
provided  to  give  an  idea  of  the  relationship. 


C  0.5  1<U 


Figure  9.  plotted  as  a  function  of  SNR^ ,  where  is  normalized  to  the  natural  frequency. 

(From  Ref.  [2].) 


A  phase  lock  loop  is  a  feedback  system  designed  to  lock  onto  the  frequency  and  phase  of 


the  input  reference  signal.  This  is  done  using  a  phase  detector,  a  loop  filter,  and  a  numerically 
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controlled  oscillator.  This  chapter  described  the  operation  of  these  components  and  derived  a 
transfer  function  for  the  system.  Using  this  transfer  function  a  typical  time  response  for  the 
system  was  plotted.  Important  figures  of  merit  were  defined  and  derived  for  the  PLL.  These 
performance  measures  are  lock  range,  lock  time,  pull-in  range,  pull-in  time  and  pull-out  range. 
Finally  a  section  describing  the  operation  of  the  PLL  in  a  noisy  environment  and  how  a  noise  can 
affect  the  figures  of  merit  of  a  system  was  included. 

The  objectives  of  this  research  are  to  implement  the  phase  lock  loop  using  fixed-point 
arithmetic.  In  order  to  do  this  a  study  of  fixed-point  arithmetic  is  needed.  The  following  chapter 
introduces  fixed-point  arithmetic  and  some  of  the  considerations  that  need  to  be  looked  at  when 
using  this  type  of  representation. 
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III.  FIXED-POINT  ARITHMETIC 


Real  world  values  are  approximated  to  nearly  any  degree  of  error  in  eomputers  using 
floating-point  arithmetie.  While  the  preeision  afforded  by  floating  point  arithmetie  is 
advantageous  in  numerous  applieations,  the  speed  at  whieh  a  eomputer  ean  ealeulate  and 
manipulate  floating-point  numbers  is  a  disadvantage.  Binary  fixed-point  arithmetie  takes 
advantage  of  the  proeessor  shift  instruetion  when  multiplying  or  dividing  by  two  to  speed  up 
arithmetie,  simply  by  representing  numbers  in  a  different  fashion.  The  seeond  advantage  to  fixed- 
point  arithmetie  is  storage  spaee.  If  the  range  of  the  real  world  values  are  known,  and  does  not 
need  32-bits  to  eover  the  range  adequately,  fixed-point  numbers  ean  be  sealed  to  eover  the  range 
and  use  perhaps  half  number  of  bits.  To  gain  an  understanding  of  how  fixed-point  arithmetie 
works,  an  introduetion  to  floating-point  representations  is  given. 


A,  FLOATING-POINT  NUMBERS 

The  IEEE  standard  754  has  dominated  most  of  today’s  proeessors  for  floating  point 
arithmetie.  It  speeifies  four  formats,  the  two  most  eommon  being  single-preeision  and  double- 
preeision.  Single  preeision  uses  32-bits  and  double  preeision  uses  64-bits.  Both  double  and 
single  preeision  formats  eontain  three  eomponents;  a  sign  bit  (S),  a  fraetion  field  (f),  and  an 
exponent  field  (e).  Eigure  10  shows  how  the  bits  are  alloeated  for  IEEE  standard  754  floating¬ 
point  format  for  both  single  and  double  preeision.  As  the  name  implies,  the  fraetion  field  is  fixed 
and  the  radix  point,  defined  by  the  exponent  field  as  deseribed  below,  is  variable.  Subsequently, 
for  small  integer  number,  most  of  the  f  field  bits  are  zeros,  henee  wasted  spaee. 

The  exponent  field  is  expressed  as  2'’“'^^  for  single  preeision  and  for  double 

preeision  where  e  is  a  variable  between  0  and  256.  The  field  seales  the  fraetion  field  and  plaees 
the  radix  point  aeeordingly.  The  range  of  numbers  able  to  be  expressed 
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SINGLE  PRECISIDV 


■'51 


(b)  DDUBLE  PREEISIOM 


Figure  10.  IEEE  754  Eloating  Point  Number  Eormat 

is  a  funetion  of  the  exponent  field,  which  ranges  from  2“'^^  -10“^*  {e  =  0)  to  2'^*  ~10^*  {e  = 
256)5. 

The  fraction  field,  as  already  alluded  to,  is  fixed,  and  therefore  is  responsible  for  the 
precision  of  the  number.  The  precision  is  defined  as  the  distance  between  two  subsequent 
representable  numbers.  This  is  simply  2“^  ,  where  f  is  the  number  of  bits  in  the  f-field.  Eor 
single  precision  floating  point,  the  precision  is  2”^^  ~  10“’ . 

B.  FIXED-POINT  NUMBERS 

As  can  probably  be  concluded,  fixed-point  numbers  fix  the  radix  point  and  the  fraction 
field  is  the  variable.  This  has  two  distinct  advantages;  1 .)  by  allowing  the  fraction  to  be  variable 
in  size,  the  storage  space  for  an  integer  can  be  tailored  to  the  range  of  values  need,  and  2.) 
arithmetic  is  essentially  integer  arithmetic  and,  when  done  in  binary,  this  amounts  to  shifts,  rather 
than  cumbersome  multiplications. 


1,  Fixed-point  Number  Representation 

Fixed-point  numbers  can  be  specialized  to  unsigned  and  signed  integers,  or  fractionals; 
however,  the  most  general  form  will  be  considered  here,  which  is  a  signed  general  fixed-point 
number,  either  an  integer,  fraction,  or  combination.  To  understand  the  usefulness  of  fixed-point 
numbers,  a  more  intuitive  look  at  decimal  fixed-point  arithmetic  will  be  presented,  followed  by 

5  A  bit  is  needed  for  the  possibility  of  exceptional  numbers  such  as  infinity  or  NaN. 
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the  extension  to  full  binary  fixed-point  arithmetie.  Fixed-point  numbers  in  a  proeessor  are 
represented  solely  as  integers.  In  order  to  properly  insert  the  deeimal  point,  a  proeessor  has  a 
syntax  for  all  numbers.  For  instanee,  if  the  proeessor  uses  a  S.3.2  syntax,  this  means  a  sign  bit 
(S),  3  integer  digits,  and  2  fraetional  digits  (realize  that  here  a  16-bit  range  is  being  used  to 
speeify  the  sign  bit  and  the  live  deeimal  digits).  The  number  238.15  in  fixed-point  S.3.2  notation 
is  represented  as  the  integer  23815.  The  original  number  is  multiplied  by  10  ,  where  the  exponent 
two  is  taken  from  the  third  field.  To  add  two  fixed-point  numbers,  eonversion  to  the  notation  of 
the  number  with  the  largest  third  field  is  required.  A  numerieal  example  will  help  to  illuminate 
the  proeess.  Suppose  we  have  the  number  above,  238.15,  in  S.3.2  notation  and  want  to  add  the 
number  145.6  in  S.4.1  notation.  The  proeessor  has  the  integer  numbers  23815  and  1456  stored, 
but  to  add  the  two,  it  must  eonvert  both  numbers  to  S.3.2.  The  whole  proeess  is  as  follows 


238.15* 

23815 

100 


100  23815  10  1456 

- = -  and  145.6*  —  = - 

100  100  10  10 

1456  23815  1456  10  38365 

10  “  100  10  10  ~  100 


whieh  is  383.65  in  S.3.2  format. 


The  payoff  eomes  when  we  extend  this  to  binary  fixed-point  numbers.  Instead  of 
multiplying  by  powers  of  10,  we  multiply  by  powers  of  2.  In  a  proeessor,  this  is  aeeomplished 
with  shifts.  Shifts  in  a  proeessor  are  extremely  fast,  thus  speeding  up  eomputationally  heavy 
applieations. 

Of  eourse  a  proeessor  does  not  keep  traek  of  the  format  as  deseribed  above,  but  the  aetions 
are  equivalent.  A  fixed-point  number  is  represented  by 

V  =  SQ  +  B  (3.1) 


where  V  is  the  real  world  number,  B  is  the  bias,  whieh  is  used  for  signed  numbers,  Q  is  the 
integer  that  eneodes  V,  and  S  =  2^  is  the  slope  or  sealing  of  the  number.  The  sealing  determines 
where  the  radix  point  belongs. 
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2.  Precision  and  Range  in  Fixed-point  Arithmetic 

The  range  of  a  number  gives  the  upper  and  lower  limits  of  representation,  while  the 
preeision  or  resolution  gives  the  distance  between  two  successive  numbers  in  the  representation. 
In  fixed-point  numbers  these  are  a  function  of  the  scaling,  S,  and  the  number  of  bits  used  or  length 
of  V.  The  range  and  precision  of  a  number  are  in  constant  competition.  For  a  large  range  the 
resolution  must  be  reduced.  Conversely,  for  a  high  resolution  the  range  is  decreased.  To 
understand  this  relationship,  calculations  of  both  of  important  parameters  will  be  made.  The 
terminology  that  is  used  will  be  a  high  or  large  resolution  is  desirable,  which  is  equivalent  to  a 
large  or  high  precision.  Often  the  term  large  precision  indicates  that  the  distance  between  two 
successive  representable  numbers  is  large.  For  the  purposes  of  this  research,  a  large  precision 
indicates  a  high  degree  of  accuracy,  and  a  low  or  small  precision  indicates  a  low  degree  of 
accuracy. 

Precision  is  determined  by  the  least  significant  bit  (LSB),  more  specifically,  what  effect 
changing  the  LSB  from  a  zero  to  a  one  has  on  the  number.  Changing  the  LSB  is  the  smallest 
increment  that  can  be  made,  thus  the  effect  of  changing  the  LSB  determines  the  precision.  So 
what  does  the  LSB  represent?  The  scaling  factor,  S,  determines  where  the  radix  point  belongs, 
hence,  determines  what  the  LSB  represents.  Suppose  E  =  0,  i.e.,  5  =  2°.  This  implies  the  radix 
point  is  to  the  right  of  the  LSB,  and  thus  a  change  from  a  zero  to  a  one  in  the  LSB  results  in  a 
change  in  the  number  of  1.  The  resolution  of  this  scaling  is  1.  For  a  better  resolution,  E  is 
required  to  be  negative.  For  E  =  -S ,  the  resolution  is  2~^  =  0.00390625 .  Generally,  the 
precision  is  defined  as  the  scaling  5  =  2^. 

This  increased  precision,  however,  comes  at  a  price.  Assume  we  are  using  n  bits  to 
represent  the  fixed-point  numbers.  This  means  that  a  total  of  2"  possible  numbers  can  be 
represented.  It  should  be  clear  that  if  a  finite  number  of  values  can  be  represented,  the  greater  the 
precision,  the  lesser  the  range.  To  determine  the  range,  simply  multiply  the  precision  by  2" 
which  gives  the  total  range  of  numbers  available.  Recall  the  bias,  B,  is  simply  a  sliding  window 
for  the  range  of  numbers.  Hardware  represents  numbers  from  0  to  2" ,  thus  for  signed  numbers  a 
negative  bias,  B  =  -2"“'  is  used,  which  slides  the  range  between  -2"“'  and  2"“' .  Table  1  shows 
the  range  and  precision  of  a  16-bit  signed  and  unsigned  fixed-point  numbers  for  various  scalings 

using  a  zero  bias  for  the  unsigned  case  and  B  =  2"“'  bias  for  the  signed. 
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Scaling 

Precision 

Range  of  Signed 
Values  (low, high) 

Range  of  Unsigned 
Values  (low, high) 

2"' 

0.125 

-4096,4095.875 

0,8191.875 

2-^ 

0.0625 

-2048,2047.9375 

0,4095.9375 

2"' 

0.03125 

-1024,1023.96875 

0,2047.96875 

2"' 

0.015625 

-512,511.984375 

0,1023.984375 

2"’ 

0.0078125 

-256,255.9921875 

0,511.9921875 

2“'* 

0.00390625 

-128,127.99609375 

0,255.99609375 

2-10 

0.00097656 

-32,31.999023 

0,63.999023 

2-12 

0.000244140 

-8,7.9997559 

0,15.9997559 

2"'' 

0.000030518 

-1,0.99996948 

0,1.99996948 

Table  1 .  Range  and  Preeision  of  a  16-bit  Fixed-point  Data  Type 


3.  Errors  in  Fixed-point  Numbers 

Due  to  the  relationship  between  range  and  preeision  in  fixed-point  arithmetie,  errors  oeeur 
between  real  world  values  and  fixed-point  representations  of  those  values.  Therefore,  it  is 
important  to  know  the  ranges  and  values  of  the  required  numbers  and  how  important  aeouraey  is 
when  using  fixed-point  arithmetie.  Errors  ean  oeeur  in  essentially  two  ways,  the  value  being 
represented  is  outside  the  range  of  the  sealing  used  and  error  due  to  a  too  low  resolution. 

The  first  kind  of  error  is  often  easier  to  find.  In  an  applieation,  an  out-of-range  error, 
while  not  neeessarily  raising  an  error  message  flag,  ean  be  found  through  analysis.  Most  fixed- 
point  proeessors,  when  eneountering  a  value  outside  the  fixed-point  range,  saturate  the  fixed-point 
representation  to  the  max  or  min  of  the  range.  This  ean  be  a  very  bad  approximation  of  the  real 
world  value.  Thus  knowing  the  range  of  values  in  an  applieation  is  vital  so  that  proper  sealing 
ean  be  used. 


29 


A  low  resolution  can  also  lead  to  errors  which  may  be  difficult  to  find.  A  processor 
representing  a  number  between  two  successive  fixed-point  numbers  does  one  of  three  things;  it 
rounds  down,  rounds  up,  or  rounds  to  the  nearest  representable  value.  To  reduce  errors  as  much 
as  possible,  rounding  to  the  nearest  representable  value  is  preferred.  This  reduces  the  error  to  half 
the  precision  and  was  the  approach  used  for  the  PLL  simulation  in  Chapter  IV.  Care  must  be 
taken  when  arithmetic  is  done  in  fixed-point  because  adding,  subtracting  and  multiplying  can 
have  dramatic  effects  on  the  size  of  numbers.  Subtracting  two  similar  numbers  can  result  in  a 
number  smaller  than  the  precision,  hence  equates  to  zero.  Equally  troubling  is  adding  a  small 
number  to  a  large  number.  The  larger  number  scaling  must  be  used  because  the  larger  range  is 
required,  but  this  means  low  precision.  If  the  small  added  number  is  smaller  than  the  precision  of 
the  larger  number,  the  addition  will  have  no  effect.  It  will  be  equivalent  to  adding  zero. 

While  increasing  the  number  of  bits  can  increase  precision  and  range  to  any  amount,  this 
is  not  a  valid  solution  for  practicality  reasons.  In  today’s  commercial  market,  32-bit  fixed-point 
processors  are  the  largest  available.  Application-specific  DSPs  are  made  with  a  larger  number  of 
bits,  but  these  are  not  readily  available  nor  cheap.  Careful  consideration  of  application  values  and 
ranges  must  therefore  be  a  priority  when  using  fixed-point  mathematics. 


Using  a  fixed-point  arithmetic  implementation  requires  knowledge  of  the  range  of  values 
and  precision  required  for  the  application.  The  range  and  precision  is  determined  by  the  scaling 
used.  For  the  fixed-point  PLL  in  the  next  chapter,  each  input  and  output  range  was  determined 
from  the  floating-point  model,  then  appropriate  scaling  used  at  each  level.  The  next  chapter 
details  the  procedure  in  developing  the  floating-point  model  and  conversion  to  a  fixed-point 
simulation. 
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IV.  FIXED-POINT  PEL  SIMULATION 


The  PLL  simulation  developed  in  this  researeh  used  Mathwork’s  Simulink  software 
paekage.  The  reasons  for  ehoosing  Simulink  for  this  applieation  are  numerous.  This  software  is 
designed  for  eontrol  systems  and  has  an  imbedded  fixed-point  bloekset.  This  feature  negated  the 
requirement  to  generate  eode  to  simulate  fixed-point  arithmetie.  Another  desirable  feature  of 
Simulink  is  the  Real  Time  Workshop,  whieh  generates  C  eode  from  a  simulation  model.  This 
feature  will  allow  easy  implementation  on  an  FPGA  as  C  ean  be  easily  eonverted  to  VHDL,  the 
language  required  for  programming  an  FPGA.  Finally,  the  Simulink  software  was  available  and 
familiar.  There  was  no  need  to  buy  expensive  modeling  or  simulation  software,  and  beeause  the 
software  was  familiar,  learning  a  whole  new  language  and  set  of  proeedures  was  negated. 

A.  PROCEDURE  FOR  MODEL 

The  proeedure  for  building  the  model  was  first  to  eonstruet  a  floating-point  simulation  in 
MATLAB  using  the  design  from  Referenee  [1].  Flaving  a  working  floating-point  model  gives 
expeeted  results  of  a  PLL  simulation  and  a  eomparison  for  the  fixed-point  model.  The  next  step 
was  to  implement  the  floating-point  model  in  Simulink,  the  bloek  diagram  based  simulation 
software  by  Mathworks.  This  done,  the  Simulink  model  was  eonverted  to  a  fixed-point  version 
where  analysis  ean  begin. 

1.  Floating-point  MATLAB  Simulation 

Any  software  simulation  that  takes  a  signal  as  input  must  rely  on  sampling  that  signal  and 
using  the  diserete  data  points.  Thus  to  the  software,  a  veetor  of  data  of  a  20  hertz  signal  looks  no 
different  than  a  veetor  of  data  bits  of  a  200  kHz  signal.  As  long  as  the  sampling  frequeney  is 
known  and  is  larger  than  the  Nyquist  limit,  the  aetual  frequeney  of  the  input  signal  ean  be 
arbitrary.  The  veetor  representing  the  input  signal  is  a  veetor  of  numbers,  regardless  of  the  aetual 
frequeney  of  the  analog  signal  it  approximates. 
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a.  Input  signal  Assumptions 

For  the  purposes  of  the  simulation  the  center  frequency,  Fg ,  of  the  NCO  was  taken 
to  be  2  kHz.  This  just  means  the  expected  frequency  of  the  input  signal  is  “around”  2  kHz.  A 
sampling  frequency  ,  of  40  kHz  was  used.  While  a  sampling  frequency  this  high  is  not 

required  for  a  2  kHz  signal,  the  resulting  output  is  much  easier  to  analyze  visually.  The  input  and 
output  amplitudes  of  the  signals  were  unity,  with  no  DC  offset.  This  assumption  is  justified  by 
realizing  that  the  input  signal  can  simply  be  normalized  by  its  amplitude  if  other  than  one,  and 
any  DC  offset  can  be  removed  by  filter  prior  to  input  into  the  PLL.  A  Hilbert  transform  is 
performed  on  the  real  input  signal  creating  a  complex  signal  with  in-phase  and  quadrature-phase 
components.  The  signal  was  initially  given  a  random  frequency  ,  within  2.5%  of  the  center 
frequency  of  the  NCO.  This  assured  the  input  signal  to  be  within  the  lock  range.  It  was  also 
given  a  random  phase  between  0  and  7t .  The  signal  to  noise  ratio  (SNR)  at  the  input  to  the  PLL 
was  assumed  to  be  SNR^  =15  dB .  This  noise  level  assured  a  lock  was  possible,  but  also  gave  an 
idea  of  how  noise  affected  the  system.  A  pre-filter  ,  of  1  kHz  around  the  expected  input  signal 
of  2  kHz  was  assumed.  This  would  leave  the  input  signal  around  2  kHz  undistorted,  but  reduce 
the  effects  of  noise.  Summarizing  the  input  signal  assumptions, 

•  Fg  =  2  kHz  . 

•  F^  =  20  kHz  . 

•  F,.  =  1  kHz  . 

•  Wj  (t)  =  cos{27r F^t  -I-  po)  -i-  i  sm{27r F^t  -i-  po)  -i-  n{t) . 
where  F^.  g  [1950,2050],  poe  [0,.7r],  and  n{t)  is  the  noise  with  SNR^  =15  dB  . 

b.  Determination  of  PLL  Parameters 

The  first  step  in  developing  an  algorithm  for  a  software  PLL  is  to  determine 
appropriate  parameters  depending  on  the  application  and  situation.  Typical  communications 
signals  have  an  SNR  of  15  dB,  thus  we  are  assured  that  the  PLL  will  lock,  although  this  will  be 
verified  once  all  parameters  are  determined.  With  the  noise  bandwidth,  not  a  factor, 
determining  the  lock  range  is  required.  If  the  expected  input  signal  has  a  frequency  of 
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around  2  kHz,  a  fair  estimate  is  that  the  actual  frequency  of  the  signal  will  not  vary  by  more  than 
5%.  This  means  we  need  a  lock  range  of  100  Hz,  or  >  200;r  rad/s .  Setting 
=  200;r  rad/s  and  using  equation  (2.28),  the  natural  frequency  in  terms  of  the  damping  factor 
is  easily  determined. 

rad/s.  (4.1) 

From  control  theory,  the  transfer  function  of  the  system  is  optimally  flat  for  ^  =  0.707 ,  but  from 
an  analysis  of  noise,  the  noise  bandwidth  as  a  function  of  ^  is  flattest  for  ^  =  0.5 ,  thus  a 
compromise  of  ^  =  0.6 .  With  this  information  the  resulting  natural  frequency  is 
=  166.66.7r  rad/s  . 


Both  PD  gain  and  loop  filter  gain  were  set  equal  to  one.  The  final  parameters  to  be 
determined  were  the  lead  and  lag  constants  of  the  loop  filter.  The  transfer  function  of  the  loop 
filter  was  given  in  Equation  (2.5).  To  determine  the  constants,  the  transfer  function  needs  to  be 
related  to  the  time  domain  input-output.  To  do  this,  the  transfer  function  of  Equation  (2.5)  is 
converted  to  its  equivalent  z-transform.  This  is  done  using  the  bilinear  transform,  such  that. 


where. 


F{z)  = 


h+biZ~' 
l  +  a,z”' 


a,  =  -1 


b 


1  +  - 


tan 


r  rj.  \ 
^2^2  j 


(4.2)6 


6  For  a  complete  discussion  on  bilinear  transforms  see  [4]. 
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and  T^=\l . 

The  inverse  of  Equation  (4.2)  must  be  found  in  order  to  derive  the  eorresponding 
difference  equation  which  will  enable  the  loop  filter  to  be  implemented  in  software.  Recall  the 
transfer  function  in  the  z-domain,  F{z) ,  is  the  input  divided  by  the  output.  From  Figure  3,  the 
input  to  the  loop  fdter  in  the  sampled  digital  world  is  (n)  and  the  output  is  Uj  (n) .  Equation 
(4.2)  can  be  written 

Uj.{z)  =  U,{z)Fiz).  (4.3) 

Substituting  in  F{z), 


Uj.  iz)  =  U,  iz)  ^  Uf  (^)(1  +  (^X^o  +  b,Z-' ) ,  (4.4) 

l  +  a,z 

and  taking  the  inverse  z-transform  back  to  the  time  domain,  the  time  difference  equation  becomes, 

Uj (n)  =  -a^Uj (n  - 1)  +  (n)  +  (n-\) .  (4.5) 

Equation  (4.5)  lends  itself  well  to  a  software  implementation. 

The  topic  of  noise  and  whether  a  sufficient  SNR  at  the  input  signal  is  present  was 
glossed  over  previously.  Because  of  the  high  noise  tolerance  of  PEL’s,  high  SNR’s  is  more  of  an 
issue  in  decoding  and  demodulating  the  signal  than  it  is  in  locking  onto  the  signal.  For  this  reason 
SNR’s  for  locking  are  exceeded  to  get  better  demodulation  and  bit  error  rate  (BER)  curves. 
However,  this  topic  will  be  addressed  now. 

Recall  that  the  SNR  at  the  input  of  the  PEL,  SNR^ ,  is  related  to  SNR^ ,  the  SNR  at 

the  output  of  the  PEE  by  Equation  (2.34).  This  relation  is  a  function  of  two  parameters:  the 
prefilter  bandwidth  and  the  noise  bandwidth.  Assuming  a  typical  prefilter  bandwidth  of  1000  Hz, 
the  only  variable  in  the  relationship  is  the  noise  bandwidth,  .  The  noise  bandwidth  is  related  to 

the  natural  frequency  and  the  damping  factor  by  Equation  (2.33).  This  equation  results  in 
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=  63.5  Hz  .  Equation  (2.34)  gives  the  output  SNR^  =  23.9  dB .  Experimentation  has  shown 
that  fee’s  will  lock  with  SNR^  >  6  dB ,  thus  the  PEE  with  the  above  parameters  should  be  well 
with  the  locking  limits  with  regards  to  noise.  The  simulation  parameters  found  are  listed  below 
for  ease  of  reference. 

•  =166. 66.7r  rad/s. 

•  ^  =  0.6. 

•  Ao)^  =  IOOtt  rad/s  (100  Hz)  . 

•  T^=  3.64e~\  =  2.3e”^ . 

•  a,  =-l,  ho  =635.1475,  h;  = -621.4397 . 


c.  The  MATLAB  Simulation 

With  the  parameters  established  above,  the  algorithm  needed  to  be  determined. 
The  feedback  loop  of  the  PEL  was  done  using  a  simple  MATLAB  for-end  loop.  This  was  done 
for  convenience  because  an  input  signal  of  a  given  length  was  assumed  for  the  simulation,  thus 
the  number  of  times  required  to  “feedback”  was  known.  This  simulated  the  fact  that  real  radio 
signal  will  be  sent  in  “blocks”  or  “packets”  in  which  the  length  of  each  packet  is  known.  Inside 
the  loop  are  the  three  blocks  of  the  PEE. 

The  first  block  is  the  phase  detector  (PD).  This  block  was  implemented  exactly  the 
way  the  PD  in  Chapter  2  was  described  by  using  the  in-phase  and  quadrature-phase  components 
of  the  signal  and  some  trigonometric  properties.  Instead  of  continuous  signals  being  multiplied 
and  manipulated,  the  sampled  data  at  corresponding  time  intervals  was  manipulated.  The  output, 
(n)  ,  of  the  PD  is  the  scalar  number  representing  the  phase  error  at  time  n. 

The  signal  (n)  was  then  sent  to  the  loop  filter  to  get  rid  of  the  higher  frequency 
terms  from  the  PD.  The  loop  filter  was  implemented  directly  from  Equation  (4.5),  producing  the 
output  signal  Uj-  (n) ,  which  is  the  error  in  the  two  signals  at  the  PD. 

The  final  block  in  the  loop  is  the  NCO.  This  block  needed  to  take  the  output  signal 
of  the  loop  filter  and  make  a  correction  to  its  oscillation  depending  on  the  magnitude  of  (n) . 
From  the  discussion  of  the  NCO  in  Chapter  2.A.3,  the  phase  of  the  continuous  time  signal  is 
given  by  Equation  (2.7).  In  the  discrete  realm,  the  phase  change  of  ^2(0  can  be  determined  by 
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(4.6) 


A02  —  {0)q  +  (n))T’^  . 

where  0)^  is  the  eenter  frequeney  of  the  NCO  in  rad/s.  With  the  ehange  in  phase  known,  the  total 
phase  of  the  signal  ean  be  determined  by 

6>2 (n  + 1)  =  6»2  (n)  +  (0)^  +  K^u ^  .  (4.7). 


If  02  is  initializd  to  02(O)  =  O5  this  computation  beeomes  possible.  Notiee  that  Equation  (4.7)  is 
iterative  and  theoretieally  02  can  continue  to  grow,  thus  a  simple  if-then  statement  is  used  to 
bound  02  between  -tt  and  tt  .  This  is  done  by  subtraeting  Itt  from  02  whenever  02>7[. 
Beeause  the  output  of  the  NCO  is  required  to  be  the  in-phase  and  quadrature-phase  eomponents,  a 
look-up  table  is  used  to  get  the  desired  output. 


U2{n)2  =  oos(02(^)) 
^2(^)2  =  sin(02(^))- 


(4.8). 


This  signal  is  now  sent  to  the  PD,  and  the  loop  is  eomplete.  A  eomplete  listing  of  the  eode  used 
for  the  MATLAB  simulation  is  given  in  Appendix  A. 

While  a  full  analysis  of  the  performanee  of  the  above  PEL  is  not  the  goal  nor 
helpful  in  this  researeh,  eonfirming  the  PEE  algorithm  operated  eorreetly  is  required.  The  input 
signal  used  had  the  assumptions  deseribed  in  Subseetion  4.A.l.a.  However,  instead  of  a  random 
frequeney,  an  input  signal  of  2100  Hz  is  used.  This  signal  represents  a  5%  error  in  the  eenter 
frequeney  of  the  model  and  the  input  signal,  whieh  is  the  worst  ease  expeeted.  Eigure  1 1  shows 
the  results  of  the  simulation.  The  output  frequeney  of  the  NCO  is  shown  in  (a)  and  the  phase 
error  of  the  two  signals  is  shown  in  (b).  The  output  ean  be  seen  to  oscillate,  but  is  damped  and 
settles  on  the  eorreet  frequeney.  The  error  plot  ean  also  be  seen  to  oseillate,  but  dies  out  and  goes 
to  zero. 

Evaluating  look-time,  reoall,  =  2.;r/<y„  =  0.0120  s,  whieh  from  the  plots  in 

Eigure  1 1  appear  to  be  aoourate.  The  oscillations  have  died  out  at  this  time  and  look  has  been 
aohieved. 
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Figure  11.  (a)  Frequency  Time  Response  of  PLL.  (b)  Phase  Error  Time  Response  of  PEL. 

Similar  evaluations  on  the  pull-in  range,  which  is  theoretically  infinite,  proved  to 
be  as  accurate  as  testable  in  the  simulation.  The  pull-in  time  of  the  system  is  exponential  with  the 
initial  frequency  displacement.  Proving  this  is  the  case  is  not  plausible  for  all  ranges  of 
frequencies.  Flowever,  given  an  initial  frequency  larger  than  the  lock  range  of  the  system,  the 
PEE  indeed  pulled  the  NCO  output  to  the  larger  frequency  offset  and  eventually  did  achieve  lock. 
As  seen  in  Eigure  12,  the  output  of  the  PD  eventually  reaches  zero  after  going  through  several 
oscillations.  This  figure  was  generated  with  an  input  signal  of  2350  Flz,  which  is  an  offset  from 
the  center  frequency  of  the  system  of  350  Flz,  clearly  outside  the  lock  range  of  100  Hz.  From 
Equation  (2.31),  the  pull-in  time  of  the  system  should  be  0.0346  seconds.  This  is  a  conservative 
time  estimate.  This  can  be  verified  from  the  plot,  where  the  oscillations  are  completely  damped 
out. 
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Figure  12.  Time  Response  of  PLL  for  an  Input  Frequency  Larger  than  the  Lock  Range,  (a). 

Output  of  System  Demonstrating  the  Signal  Locks  on  a  Frequency  of  2350  FIz.  (b).  Plot  of  the 

Error  of  the  System. 

This  concluded  the  floating  point  MATLAB  simulation.  The  next  step  was  to 
build  this  system  in  the  control  simulation  language  of  Simulink. 

2,  Fixed-Point  Simulink  Model 

The  fixed-point  blockset  in  MATLAB  is  a  toolbox  used  in  Simulink.  Simulink  is  a 
supplement  to  MATLAB  for  building  models  and  running  simulations,  thus  ideal  for  the  purposes 
of  the  PLL.  It  uses  blocks  as  functions  and  mathematical  operators.  Implementing  the  functional 
blocks  of  the  PLL  was  a  matter  of  translating  the  MATLAB  code  to  blocks  in  Simulink.  The 
model  can  be  seen  in  Appendix  A. 

The  Simulink  model  was  built  using  the  same  parameters  as  the  MATLAB  simulation.  As 
expected,  its  performance  was  identical  to  the  performance  of  the  MATLAB  coded  model,  thus 
no  results  will  be  provided  here. 
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3.  Fixed-Point  Model 


The  point  of  this  thesis  was  to  establish  the  plausibility  of  a  fixed-point  PLL  in  software 
and  to  determine  what  type  of  performanee  ean  be  expected  from  such  a  system.  This  was  to  be 
achieved  by  building  a  model  and  running  simulations  to  analyze  the  results.  In  building  the 
fixed-point  model,  the  same  PLL  parameters  as  the  floating-point  model  were  used. 

As  will  be  the  case  throughout  the  fixed-point  process,  the  design  problem  faced  at  the 
input  and  output  of  each  block  is  determining  an  appropriate  scaling  for  the  fixed-point  number  so 
that  a  large  enough  range  is  used,  without  introducing  unacceptable  inaccuracies  due  to  a  lower 
precision.  As  will  be  explained  in  Chapter  4.C,  a  16-bit  fixed-point  implementation  was  used. 
Essentially,  the  reason  for  this  was  while  an  11 -bit  representation  was  sufficient  to  run  correctly, 
typical  FPGA’s  and  DSP  chips  use  8-bit,  16-bit,  or  32-bit  fixed-point  arithmetic.  The  8-bit 
representation  didn’t  give  enough  precision  for  the  system  to  lock  with  any  reliability,  thus  16-bits 
were  used. 

The  procedure  to  convert  the  floating-point  Simulink  model  to  a  fixed-point  model 
involved  looking  at  each  block  and  determining  the  range  of  the  signal  at  the  input  of  the  block 
and  the  range  of  the  signal  at  the  output  of  the  block.  With  this  information,  a  scaling  appropriate 
for  the  range  can  be  determined  and,  subsequently,  the  precision  can  be  found  at  that  scaling.  If 
the  precision  was  too  large,  an  alternate  method  must  be  used,  or  an  increase  in  the  number  of  bits 
used. 

To  begin,  the  input  signal  to  the  PLL  needed  to  be  converted  to  its  fixed-point  equivalent. 
Recall  the  input  signal  was  normalized  to  ±1.  From  Table  1,  the  best  fit  for  this  range  of  numbers 
is  using  a  scaling  of  2“'^ .  Because  the  range  at  this  scaling  is  exactly  the  expected  range  of  the 
input,  any  noise  could  be  amplified,  hence  a  slightly  larger  range  was  used,  specifically  a  scaling 
of  2“'"^.  This  gives  a  precision  of  6.1035e“%  which  is  still  below  the  noise  of  the  input  signal, 
thus  is  acceptable  resolution. 

With  the  input  signal  converted  to  fixed-point,  the  logical  next  step  was  to  convert  the  PD 

to  a  fixed-point  equivalent.  The  block  diagram  of  the  PD  is  shown  in  the  Appendix  in  Figure  A2 

(duplicated  on  the  next  page  for  convenience).  Tracing  the  connections,  one  can  see  that  this  is 

equivalent  to  the  MATLAB  code  implementation  in  Appendix  A.  The  PD  takes  the  input  signal 
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and  multiplies  it  by  the  in-phase  and  quadrature -phase  components  of  the  NCO.  In  the  block 
diagram  of  Figure  A2,  In2  (short  for  input  2)  is  the  phase,  >  which  goes  into  a  look-up  table 
for  sine  and  cosine  producing  the  in-phase  and  quadrature-phase  signals.  The  multiplication 
mentioned  above  multiplies  two  signals  with  a  range  of  values  between  -1  and  +1,  thus  the  range 
of  the  output  of  this  block  must  be  within  the  range  of-1  and  +1.  As  can  be  seen  in  Figure  A2, 
below  the  multiplication  blocks  Productl  through  Productd,  one  can  see  the  scaling  is  again  2“'"^ , 
reflecting  the  range  required  at  the  output  of  this  block. 


Ftm 


Figure  A2.  Phase  Detector  Simulink  Model.  (Duplicated  From  the  Appendix  for  Ease  of 

Reference.) 

To  belabor  the  point,  the  next  blocks  in  the  PD,  Suml  and  Sum2,  simply  add  or  subtract  the 
output  of  the  multiplication  blocks,  producing  a  range  of  -2  to  +2,  which  is  the  range  of  the 
scaling  2“'"^ .  As  can  be  surmised,  this  process  of  determining  the  range  at  every  block  needs  to  be 
done,  and  appropriate  scaling  of  the  output  determined.  Care  must  be  taken  in  ensuring  that  when 
applying  the  best  possible  range  to  the  output,  a  resolution  problem  isn’t  introduced.  This  issue 
will  be  discussed  in  Section  4.C.  and  is  the  reason  an  8-bit  implementation  was  not  possible.  The 
full  block  diagram  of  the  Simulink  PEL  is  included  in  Appendix  A,  with  scaling  at  each  block 
indicated. 
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The  LF  was  designed  using  the  1/z  bloek.  Referring  back  to  the  z-transform  of  the 
difference  equation  established  for  the  MATLAB  simulation  and  rearranging  equation  (4.4), 

Uf{z)  =  -a,z~'U ^  (z)  +  +  b,z~'U ^  (z) ,  (4.9) 

where  is  the  unit  delay  operator.  This  equation  is  the  one  used  in  the  Simulink  model  for  the 
loop  filter. 

The  NCO  was  implemented  similarly  using  constant  blocks  and  multiplier  blocks, 
ensuring  the  range  of  the  inputs  and  outputs  were  scaled  properly.  The  complicated  circuitry  in 
the  NCO  is  due  to  the  fact  that  0^  grows  each  increment  and  must  be  reduced  to  -7t  to  7t .  This 
was  done  using  a  comparator  block  and  a  switch.  The  output  of  the  switch  is  the  top  input,  which 
is  simply  9^  for  ^2  <  •  If  ^2  >  ’  the  comparator  resulted  in  zero,  the  switch  outputs  its  third 

input  vice  its  first,  which  is  9^-ln  . 

B,  PERFORMANCE  AND  ANALYSIS  OF  SIMULINK  PLL  MODEL 

The  performance  of  the  model  was  consistent  with  the  MATLAB  simulation  considered  in 
Subsection  4.A.I.C.  For  high  SNR  values,  the  fixed-point  model  performed  identically  to  the 
floating-point  model  in  lock  range,  lock  time  and  pull-in  time.  The  error  of  the  output  signal  due 
to  the  quantization  of  the  fixed-point  signal  became  evident  in  the  fixed-point  model.  This  error  in 
effect  reduces  the  SNR  of  the  output  signal. 

1.  Performance  Measures  of  the  Fixed-point  Model 

Figure  13  shows  the  time  response  of  the  output  of  the  PD,  corresponding  to  the  error  in 
the  PLL.  It  is  clear  that  the  error,  6^ ,  dies  to  zero,  thus  confirming  the  PLL  achieves  lock.  To 
generate  this  plot  a  reference  frequency  of  2100  Hz  was  used. 
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Figure  13.  (a).  Plot  of  Fixed-point  Model  NCO  Output,  (b).  Plot  of  Fixed-point  Model  Error. 

To  establish  the  performance  of  the  PEL,  the  theoretical  values  of  the  performance 
measures  will  be  computed  and  then  compared  to  the  results  obtained  by  the  simulation.  Erom 
Section  2.C,  five  key  performance  measures  were  derived.  These  are  lock-range,  lock-time,  pull- 
in  range,  pull-in  time,  and  pull-out  range.  Eor  the  parameters  used  in  the  simulation,  the 
following  theoretical  values  for  the  fixed-point  PEE  are  calculated. 

1.  =  200.;r  rad/s . 

2.  =0.0120  s. 

3.  AcOp  =  oo . 


4. 


AO)q 


5.  AcOpQ  =1508  rad/s . 


Eor  the  testing  of  the  fixed-point  PEE,  a  worst-case  scenario  input  was  adopted.  With  this 
in  mind,  the  random  phase  was  dropped  in  favor  of  an  artificially  placed  phase  offset  of  90 
degrees.  This  represents  the  worst  possible  phase  offset  and,  thus,  the  most  conservative 
measures  of  performance. 
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a.  Lock  Range 

The  lock  range  of  the  PLL  model  is  the  range  of  frequencies  at  which  the  system 
locks  in  a  single  operating  cycle.  This  corresponds  to  one  period  of  the  natural  frequency  .  To 

determine  this  value  experimentally,  the  input  signal  frequency  was  iteratively  increased  until  the 
system  did  not  lock  in  a  single  cycle.  This  established  a  bracket,  or  upper  and  lower  bound  on  the 
lock  range.  By  successively  narrowing  this  bracket  down,  the  lock  range  was  determined.  Using 
this  procedure  the  lock  range  was  found,  Ao)^  =5\%7r  rad/s,  a  value  over  twice  as  large  as  the 
theoretical  value  of  200.^  rad/s . 

While  a  large  lock  range  is  desirable,  some  explanation  as  to  why  it  is  larger  than 
predicted  is  required.  There  are  two  reasons  this  is  the  case.  First,  when  determining  the  lock 
range,  an  approximation  on  the  gain  of  the  loop  filter  at  the  frequency  offset  was  used.  The 
assumption  was  that  the  lock  range  was  larger  than  the  comer  frequencies  of  the  two  time 
constants,  and  thus  the  comer  frequency  was  used  as  an  approximation.  The  lock  range  is  indeed 
larger  than  the  corner  frequencies,  but  with  the  non-linearity  of  the  loop  filter  transfer  function,  a 
tighter  bound  is  difficult  to  calculate,  thus  the  overly  conservative  loop  filter  comer  frequency  is 
used. 

Second,  the  fact  that  noise  was  not  a  factor  in  the  design  of  this  PLL,  the  noise 
bandwidth  is  very  small,  which  means  the  noise  is  not  hampering  the  locking  of  the  device.  This 
results  in  a  larger  lock  range  than  theoretically  calculated. 

b.  Lock  Time 

The  lock  time  can  be  visually  seen  in  the  plot  of  Figure  13,  but  for  a  more  accurate 
determination,  the  discrete  output  values  were  analyzed  and  the  lock-time  determined  to  be  when 
the  output  entered  and  remained  within  2%  of  the  final  value  of  the  system.  This  is  an  accepted 
percentage  for  the  settling  time  of  a  system  in  control  theory.  This  is  equivalent  to  the  time 
required  for  the  error  to  die  down  and  remain  below  0.02.  However,  the  actual  analysis  is  more 
involved.  With  the  introduction  of  noise,  a  large  noise  value  at  any  given  sample  can  throw  the 
output  outside  the  2%  range,  even  though  the  output  due  to  the  system  is  within  the  2%  range. 
Examining  the  numerical  values  at  each  sample  and  finding  where  the  output  stayed  below  0.02, 
the  lock-time  was  found  to  be  at  0.0128  s.  Comparing  this  value  to  the  theoretical  value  of  0.0120 
s,  the  fixed-point  representation  does  not  degrade  performance  in  this  area  significantly.  It  must 
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be  understood  also  that  the  theoretieal  value  of  0.0120  is  assuming  no  noise,  and  that  any 
significant  spike  in  the  noise  at  the  output  of  the  system  could  effect  the  settling  time  of  the 
system  when  analyzing  the  output  at  each  time  sample.  However,  it  should  be  noted  that  running 
the  simulation  several  times,  the  above  settling  time  of  0.0128  was  the  average  of  all  simulations. 

c.  Pull-in  Range 

The  pull-in  time  for  the  second-order  PLL  is  theoretically  infinite,  but  an  analysis 
of  the  pull-in  range  of  the  fixed-point  system  needs  to  be  done.  At  some  point,  the  frequency 
offset  will  get  so  large  that  the  fixed-point  precision  or  range  will  be  exceeded  and  critical  errors 
will  be  introduced.  Incrementally  increasing  the  input  frequency,  the  fixed-point  PLL  eventually 
obtained  lock  for  a  max  frequency  of  2830  Hz.  This  value  represents  a  41.5%  error  from  the 
center  frequency  of  the  PLL. 

The  limitation  seems  to  come  from  the  loop  filter.  The  range  or  precision, 
depending  on  the  scaling  used,  is  exceeded  at  the  output  of  the  loop  filter,  causing  the  input  to  the 
NCO  to  be  inaccurate.  The  same  limitation  was  discovered  when  attempting  to  develop  an  8-bit 
fixed-point  PLL  and  will  be  discussed  in  more  detail  in  Section  4.D. 

d.  Pull-in  Time 

The  pull-in  time,  Tp ,  of  the  system  is  dependant  on  the  frequency  offset  from  the 
center  frequency.  Realize  that  if  the  frequency  offset  is  inside  the  lock  range,  the  lock  time  is 
used  and  is  the  same  for  any  frequency  offset,  as  long  as  it  is  smaller  than  the  lock  range.  The 
pull-in  time  varies  with  the  magnitude  of  the  offset.  To  get  an  idea  of  the  performance  as  related 
to  pull-in  time,  five  input  frequencies  were  used  and  compared  to  the  theoretical  values  for  those 
frequency  offsets.  The  settling  time  again  was  calculated  using  a  2%  tolerance.  Table  2  shows 
the  results  for  the  five  frequencies.  The  pull-in  time  is  given  in  seconds  for  the  theoretical  as  well 
as  the  experimental  and  the  frequency  offset  is  given  in  Hz.  An  initial  offset  of  300  Hz  was 
chosen  due  to  the  fact  that  from  the  analysis  of  lock-range,  =  259  Hz  and,  thus,  a  larger 
value  than  this  was  needed.  The  pull-in  range  of  the  system  was  determined  to  be  830  Hz,  thus 
only  values  up  to  800  Hz  were  included. 


Aw  [Hz] 

Tp  [s]  (theoretical) 

Tp  [s] 

(experimental) 
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0.0254 

0.0452 

500 

0.0707 

0.0446 

600 

0.1018 

700 

0.1385 

0.0836 

800 

0.1810 

0.1114 

Table  2.  Pull-in  Times  for  Various  Frequeney  Offsets 


The  disagreement  between  the  theoretical  and  experimental  values,  especially  as 
the  frequency  offset  increases,  cannot  entirely  be  explained  by  the  approximations  used  in 
deriving  Tp  .  Another  explanation  lies  in  the  discrete  quantization  of  the  fixed-point  PLL.  The 
output  of  the  PD,  or  the  error  of  the  system,  has  a  certain  precision.  The  range  of  numbers 
required  at  that  block,  as  described  previously,  determines  this  precision.  If  the  value  at  the 
output  of  the  PD  is  smaller  than  this  precision,  the  value  is  rounded  to  zero.  By  the  nature  of  the 
second-order  PLL,  any  time  two  successive  outputs  fall  within  the  precision  limits,  the  output 
remains  at  zero.  The  above  theoretical  values  assume  an  analog  PLL,  where  precision  is  infinite. 
A  fixed-point  software  PLL  has  quantization  error,  which  in  this  case,  output  a  zero  even  though 
the  real-world  value  is  larger  than  that. 

e.  Pull-out  Range 

The  pull-out  range  is  an  important  figure  because  it  allows  the  user  to  determine 
how  large  a  frequency  step  can  be  supplied  to  the  input  without  unlocking  the  system.  This 
assumes  the  system  is  already  in  lock  and  a  frequency  step  occurs,  as  in  an  FM  system.  This 
value  was  tested  by  giving  the  PLL  a  signal,  allowing  lock  to  occur,  then  stepping  up  the 
frequency  until  it  unlocked.  Simulations  of  the  fixed-point  PLL  agreed  with  the  expected  figure, 
which  was  obtained  through  computer  modelling  in  Reference  [5].  The  fixed-point  PLL 
experimentally  had  a  pull-out  range  of  1501  rad/s  or  239  Hz,  compared  to  the  theoretical  value  of 
1508  rad/s,  or  240  Hz.  Recall  that  the  lock-range  of  the  system  was  259  Hz,  thus  if  a  frequency 
step  larger  than  the  pull-out  range,  but  less  than  the  lock  range,  the  system  will  again  lock  in  the 
lock- time,  or  0.0120  seconds  for  this  system.  If  the  frequency  step  is  larger  than  the  lock-range, 
the  system  will  acquire  lock  using  the  longer  pull-in  time  as  described  above. 

/.  Frequency  drift 
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The  magnitude  of  a  frequeney  drift  is  of  interest  because,  for  a  software  radio 
implementation  and  downloadable  software  radio  packages,  satellite  communication  is  required. 
As  neither  satellites  nor  users  are  stationary,  a  frequency  drift  can  occur  due  to  the  Doppler  shift 
in  frequency.  The  Doppler  shift  is  governed  by  two  parameters:  the  velocity  and  direction  of 
travel  relative  to  the  antenna,  and  the  wavelength  of  the  signal.  Specifically,  the  Doppler  shift  in 
frequency  is, 

f,=jcose.  (4.10) 

where  v  =  velocity,  A  =  wavelength  or  d  f^,  c=speed  of  light ,  and  0  =  direction  of  travel  from 
the  antenna.  For  a  frequency  drift  to  occur,  the  satellite  would  need  to  either  be  changing  velocity 
or  changing  direction.  With  an  orbital  satellite,  the  direction  of  movement  is  constantly  changing. 
Therefore,  an  idea  of  what  kind  of  drift  the  PLL  can  handle  is  an  interesting  figure  of  merit. 

From  Subsection  2.B.1.C,  the  rate  of  change  in  a  frequency  ramp  that  will  cause  a 
system  to  unlock  is  governed  by  Equation  (2.24).  To  test  the  fixed-point  PLL  model,  a  frequency 
ramp  was  applied  to  the  input  after  the  system  was  allowed  to  lock.  The  slope  of  this  ramp  was 
increased  in  each  run  until  the  system  was  observed  to  unlock  and  to  be  unable  to  track  the 
frequency  change.  The  rate  at  which  the  PLL  was  unable  to  track  the  frequency  ramp  was  much 
smaller  than  that  predicted  by  Equation  (2.24),  but  this  value  was  assuming  no  noise  and  an 
analog  system.  Because  of  the  non-linearity  of  a  frequency  ramp,  noise  and  quantization  error  are 
magnified.  The  results  of  the  fixed-point  PLL  are  =15,000,  compared  to  the  theoretical 

value  of  130,000.  The  PLL  could  maintain  lock  up  to  this  rate  of  change  for  about  0.25  seconds, 
corresponding  to  about  30  cycles  of  the  system.  Lrom  0.25  seconds  on,  the  large  frequencies 
were  out  of  range  of  the  fixed-point  representation.  The  PLL  maintains  lock  longer  for  a  more 
gradual  frequency  ramp  slope  (see  example  below).  The  limiting  factor  is  simply  the  range  of 
frequency  that  must  be  covered  for  the  case  of  a  frequency  drift.  Lor  a  fixed-point  PLL  with  a 
center  frequency  of  2000  Hz,  the  system  unlocked  when  the  reference  frequency  reached  around 
3300  Hz.  Lor  smaller  sloped  frequency  ramps,  the  system  remained  locked  longer,  but  once  it 
reached  the  upper  bound  of  3300  Hz,  the  system  fell  apart. 

To  understand  if  this  is  a  limitation  for  a  real  world  application,  an  example  of  a 
moving  user  is  given.  There  are  two  ways  a  frequency  ramp  can  be  seen  at  the  receiver  of  a  PLL 
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due  to  the  movement  of  a  user.  The  first  is  if  the  user  is  accelerating  in  speed.  The  second  is  if 
the  user  is  changing  direction.  The  second  case  will  be  considered  first.  The  largest  change  in 
frequency  due  to  a  change  in  direction  occurs  when  0  goes  from  1  to  -1.  This  represents  a  user 
turning  around  and  moving  in  the  opposite  direction.  The  user  will  be  assumed  to  travel  at  a 
constant  velocity  of  60  mph  (26.82  m/s).  If  the  user  is  assumed  to  turn-around  on  a  dime,  the 
result  is  a  frequency  step.  The  performance  of  a  frequency  step  was  covered  in  Subsection 
d.B.l.e,  therefore  a  circular  path  of  the  user  will  be  assumed,  which  represents  a  constant  change 
in  angle,  thus  a  frequency  ramp  as  the  input  signal  of  the  PLL.  As  shown  in  Figure  14,  the  user 
makes  a  tight  circle  with  a  radius  of  10  m.  To  get  anything  but  negligible  values,  a  high 
frequency  signal  of  100  MHz  is  used.  At  26.82  m/s  the  turnaround  would  take  them  2.34 
seconds.  The  value  0  in  Equation  (4.10)  goes  from  0  to  .tt  in  2.34  seconds.  This  represents  a 
slope  of  the  frequency  ramp  of  48  radians.  A  slope  of  48  radians  would  take  170  seconds,  or 
almost  three  minutes  for  the  PLL  to  unlock.  The  frequency  ramp  is  only  applied  for  2.34  seconds. 
The  system  can  easily  handle  a  frequency  ramp  of  this  slope,  even  with  an  elevated  frequency  of 
100  MHz. 


Ligure  14.  Vehicle  Driving  60  mph  Around  a  Tight  Curve  to  Illustrate  how  the  Doppler  Shift  can 

Cause  a  Lrequency  Ramp  of  a  Signal. 

The  second  case  is  when  the  user  is  accelerating.  Suppose  the  user  goes  from  a 
dead  stop  to  60  mph  in  3.1  seconds  (the  fastest  production  car  ever  built  is  the  Porsche  Carrera 
twin  turbo  which  was  tested  at  0-60  mph  in  3.1  seconds).  The  resulting  frequency  ramp  has  a 
slope  of  18.12  radians.  Lor  a  slope  of  A<i>  =  18.12  rad/s  ,  the  PLL  would  take  450.78  seconds  to 
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unlock.  Even  in  an  accelerating  jet,  the  frequency  ramp  caused  by  acceleration  is  well  within  the 
limits  of  the  PEL 

Eor  normal  operating  environments,  the  PEE  can  be  expected  to  stay  locked  for  a 
real-world  frequency  ramp.  If  an  environment  where  the  frequency  ramp  is  too  large,  or  applied 
for  too  long  is  encountered,  a  second  PEE  with  an  appropriate  center  frequency  could  be  used. 

2,  Analysis  of  PLL 

The  performance  measures  of  the  PEE  discussed  in  the  previous  section  give  a  good  idea 
of  the  capabilities  of  the  fixed-point  PEE,  but  two  further  items  concerning  the  PEE  need  to  be 
analyzed;  the  error  at  the  output  due  to  the  fixed-point  representation,  and  why  an  8-bit  system 
failed.  The  first  will  be  examined  first  as  it  will  shed  some  light  on  the  second. 

a.  Errors  in  the  Output  of  the  PLL 

To  get  an  idea  of  the  errors  expected  at  the  output  of  the  PEE,  an  analysis  of  the 
errors  that  could  possibly  accumulate  at  each  stage  of  the  PEE  is  required.  A  full  discussion  of 
error  analysis  and  propagation  of  errors  is  beyond  the  scope  of  this  research;  however,  to 
determine  the  size  of  errors  at  the  output  of  the  PEE,  some  error  analysis  is  required.  This  will  be 
done  by  looking  at  the  precision  at  each  stage  of  the  PEE  and  from  that  precision  a  worst-case 
error  can  be  determined.  By  finding  the  largest  possible  error  in  magnitude  at  each  block,  the 
error  of  the  sinusoidal  output  of  the  NCO  is  determined.  Errors  at  the  output  of  each  block  are 
due  to  two  factors.  The  first  is  the  error  at  the  input  of  the  block.  The  second  is  the  error  due  to 
the  fixed-point  representation  and  arithmetic  operation  of  the  block.  Adding  these  two  error 
results  in  the  error  at  the  output  of  the  block. 

The  input  signal  was  assumed  to  be  normalized  to  amplitude  of  ±1 ,  thus  the  best 
scaling  of  the  fixed-point  number  that  includes  that  range  is  5  =  2“'"^  (see  Chapter  3  on  scaling  of 
fixed-point  numbers).  This  results  in  a  precision  of  h.lOE-S  .  Recall  that  precision  is  defined  as 
the  distance  between  two  successive  representable  numbers.  If  a  “round  to  nearest” 
implementation  is  used,  the  largest  possible  error  that  can  occur  is  half  of  this  precision  or 
3.05E-5  .  The  convention  of  =  accumulated  error  will  be  adopted  for  the  remainder  of  this 
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section.  With  this  notation,  when  the  real  world  signal  is  converted  to  a  fixed-point  value, 
e^=3.05E-5. 

For  the  PD  component  of  the  PLL,  the  input  signal  is  multiplied  by  the  output  of 
the  NCO.  The  initial  value  of  the  NCO  is  zero,  but  it  is  sinusoidal  in  nature.  The  assumption  will 
be  made  here  that  the  initial  error  in  the  output  of  the  NCO  is  due  only  to  the  precision  of  the 
fixed-point  number  representation.  The  NCO  sinusoid  output  also  has  amplitude  ±1,  thus  uses 
the  same  range  and  precision.  For  multiplication,  the  error  in  the  product  is  the  sum  of  the  error 
in  the  two  multiplicands.  The  error  of  both  multiplicands  is  the  same,  thus  after  the  multiplication 
of  the  reference  signal  and  the  NCO,  the  error  due  to  the  input  becomes 
2*3.05£’-5  =  6.1035£’-5.  Adding  the  error  due  to  the  fixed-point  operation  yields  a  value  of 
e^=9.\552E-5. 

The  next  block  adds  the  in-phase  and  quadrature -phase  components  in  the  PD. 
Again,  finding  the  error  due  to  the  inputs  and  adding  this  value  to  the  error  due  to  the  fixed-point 
operation  determines  the  error.  For  addition,  the  error  due  to  the  inputs  is  the  sum  of  the  errors. 
This  results  in  a  value  of  1.831  OF' -4  .  Adding  the  resolution  error,  the  output  of  the  summation 
has  a  possible  error  of  magnitude  =  1.8310F-4-I-3.0517F-5  =  2.1362F-4  . 

The  next  step  of  the  PD  was  to  divide  the  in-phase  and  quadrature-phase  signals. 
To  determine  the  error  at  the  output  of  this  operation,  assume  first  that  the  division  operation  is 
perfect  and  that  the  output  error  is  caused  only  by  the  error  at  the  input.  This  being  the  case,  the 
error  propagates  the  same  as  multiplication.  The  resulting  error  is  the  sum  of  the  errors  at  the 
input.  Therefore  the  error  at  the  output  due  to  error  at  the  input  is  2  *  2. 1362F  -  4  =  4.2724F  -  4 . 
Unfortunately  the  errors  of  this  operation  do  not  stop  here.  If  the  denominator  in  this  operation  is 
small  compared  to  the  numerator,  a  large  number  will  result.  In  fixed-point  arithmetic,  this  means 
a  larger  range  of  numbers  needs  to  be  represented.  For  the  division  in  the  PD,  the  range  of 
numbers  was  determined  to  be  between  ±128  .  To  get  this  type  of  range,  a  scaling  of  5  =  is 
needed,  which  results  in  a  precision  of  3.906F-3  .  The  worst-case  error  at  this  stage  is  the  error 
due  to  the  input  plus  the  error  due  to  fixed-point  scaling.  Adding  these  two  errors  together, 
=2.3804F-3. 
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The  next  bloek  of  the  PD  is  the  tan  '  (x)  funetion.  This  funetion  is  asymptotieally 

bounded  by  ±n H,  thus  a  sealing  of  5  =  2“'"^ is  used.  Beeause  this  is  a  non-linear  funetion,  errors 
do  not  work  exaetly  the  same.  The  error  will  be  maximized  where  the  funetion  has  a  maximum 
derivative.  This  is  at  x  =  0  .  To  find  the  error  at  the  output,  the  funetion  must  be  applied  to  the 
error  at  x  =  0  .  More  speeifieally,  if  the  error  at  the  input  of  this  stage  is  =2.3804£'-3,  the 

represented  value  of  the  fixed-point  system  could  be  0  or  2.3804£’-3  .  If  the  tan”'  function  in 
the  system  was  perfectly  continuous,  the  error  at  the  output  could  be  as  large  as 

=tan”'(0)-tan”'(2.3804£-3) 

=  2.3804£-3. 

Including  the  error  due  to  the  fixed-point  representation  of  tan”'  function,  the  error  is 
=  2.3804£'-3-l-3.0517£'-5  =  2.4109£'-3  .  As  this  is  the  last  block  of  the  PD,  the  total  error 
of  the  phase  detector  component  is  e^  pj^  =  2.4109£’  -  3 . 

The  loop  filter  component  is  the  biggest  contributor  to  error  in  the  system.  This  is 
because  the  loop  fdter  lead-lag  constants  are  large,  thus  a  larger  range  is  needed,  reducing 
precision.  The  input  to  the  loop  filter  is  the  output  of  the  tan”'  function,  which  has  a  range  of 
±71 H .  This  signal  is  multiplied  by  the  constants  ho  and  hj .  For  the  PLL  model,  these  two  values 
were  635. 14±  0.015625  and  -621.43  ±0.015625,  respectively.  The  output  of  this  multiplication 
needed  to  include  the  entire  range  of  ±(.;r/ 2)(635.14)  ~  ±997  .  This  required  a  scaling  5  =  2”% 
which  has  a  precision  of  0.03125  or  output  error  of  0.015625 .  Adding  the  error  due  to  the  inputs 
of  the  multiplication  to  the  error  from  fixed-point  multiplication,  the  accumulated  error  becomes 
=  3. 366 IF' -2  on  each  multiplication. 

The  error  gets  worse  when  the  addition  in  Equation  (4.5)  of  the  loop  filter  is 
applied.  This  addition  has  three  terms  in  it.  The  first  term  in  Equation  (4.5)  has  the  same  range  as 
the  output  of  this  addition.  Because  of  the  feedback  of  the  system,  a  range  of  ±5000  is  needed. 
This  results  in  a  precision  of  0.25  or  error  of  0.125.  The  error  in  the  three  terms  are  added 
together  then  added  to  the  precision  of  the  output.  The  total  error  of  the  loop  filter  then  becomes, 

=2*3.3661F-2  +  0.125  +  0.125  =  0.2837.  (4.11) 
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The  first  block  in  the  NCO  performs  the  operation 

0Xn)  =  O)QTs+Uf{n)Ts.  (4.12) 

The  first  term  in  this  sum  is  a  constant.  Evaluated  to  be  0.62831,  it  is  the  center  frequency  of  the 
PEL  multiplied  by  the  inverse  of  the  sampling  frequency.  Eor  this  PEE  model 
=1/Ej  =5E-5  .  The  second  term  varies  with  the  signal  Uj  {n);  hence  a  range  of  numbers  is 

required  vice  the  constant  first  term.  The  range  of  Uf{n)  is  the  range  of  the  output  of  the  loop 
fdter,  which  was  previously  determined  to  be  ±5000 .  Thus  the  range  of  the  second  term  in  the 
sum  is  ±5000 =±0.25  .  Eor  this  small  range,  a  scaling  of  S  =  2“'^  is  used,  which  increases 
precision.  To  determine  the  error  at  the  output  of  this  multiplication,  the  error  due  to  the  precision 
of  the  output  is  added  to  the  largest  possible  error  at  the  input,  multiplied  by  .  Thus  the 

error  at  the  output  becomes, 

+i2-'-’=2.9442£-5. 

Adding  the  first  term  in  the  sum  of  Equation  (4.12),  a  new  range  is  needed  due  to  the  larger  value 
of  0.62831 .  This  larger  range  results  in  a  scaling  of  5  =  .  The  error  due  to  this  precision  is 

3.0517£'-5.  Adding  the  error  in  the  fixed-point  representation  of  0.62831,  the  error  in  the 
second  term  of  2.9442£'-5  and  the  error  due  to  the  fixed-point  operation  yields  an  error  of 
e,  =9.0476E-5. 

The  final  possible  introduction  of  errors  in  the  system  is  due  to  feedback  property 
of  the  NCO  and  the  continual  scaling  of  the  output  O^in) io  keep  it  between  -7t  and.7r.  The 
output  phase  of  the  NCO  equals,  9^{n)  =  9j^{n-\)  +  C0QT^+Uf{n-\)T^.  If  is  assumed  to 

remain  between  -7t  and  7t ,  92{n-\)  must  be  in  this  range  as  well.  The  appropriate  scaling  for 

this  sum  is  5  =  2“'^ .  This  results  in  an  error  of  6. 1035E'  -  5  .  The  last  two  terms  of  the  sum  have 
a  total  error  of  9.0476£'-5 .  Adding  the  errors  due  to  fixed-point  scaling  and  errors  due  to  the 
feedback  input,  the  resulting  error  of  ^2(^)  is  equal  to 
=9.04756E-5  +  2*6.1035E-5  =  2.1254E-4. 
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To  keep  02in)  bounded,  the  algorithm  adopted  was  to  subtraet  In:  whenever  it  got 
larger  than  n .  To  get  an  aecurate  error  of  the  system,  this  subtraetion  must  be  taken  into  aecount, 
despite  the  faet  that  it  oeeurs  infrequently.  To  subtraet  2n  from  the  phase  whenever  it  grows 
larger  than  n  requires  eonverting  2n  to  a  fixed-point  value.  The  error  in  this  eonversion  is 
6.1 03 5 £ -5.  Again  adding  the  error  due  to  the  inputs  to  the  error  due  to  the  fixed-point 
operation,  the  total  error  of  the  NCO  is  determined  to  be, 

3  •3462£-4.  (4.13). 

The  error  determined  thus  far  is  the  error  in  the  phase  of  the  NCO.  The  output  of 
the  NCO  is  a  sinusoid  (aetually  two  sinusoids,  the  in-phase  and  quadrature-phase  eomponents  of 
the  signal,  but  the  error  will  be  the  same  for  both)  with  amplitude  one.  The  sealing  used  for  this 
range  is  5  =  2“'"^ .  To  determine  the  output  error  of  this  sinusoid,  the  error  at  the  input  must  be 
applied  to  the  sinusoid  funetion  where  its  derivative  is  maximum.  For  the  sine  funetion,  this  is  at 
zero,  thus  the  error  at  the  output  due  to  the  input  error  is  sm(3.3462£'-4)-sm(0)  =  3.3462£'-4  . 

This  value  must  be  added  to  the  error  of  the  sinusoid  due  to  its  preeision.  With  a  sealing  of  2“'"^  a 
fixed-point  quantization  error  of  3. 051 VF'- 5  is  found.  Adding  these  two  errors  together,  the 
total  error  of  the  system  is  determined  to  be 

C,TOr,i=3.3462£-4  +  3.0517£-5  =  3.6513£-4.  (4.14) 

The  subseript  (1)  refers  to  the  error  after  the  first  iteration  of  the  feedbaek  loop.  At 
the  beginning  of  this  analysis,  the  NCO  error  was  multiplied  by  the  referenee  signal  input. 
Beeause  this  is  a  feedbaek  loop  and  a  start  point  had  to  be  determined,  the  NCO  input  error  was 
assumed  to  be  zero.  Now  that  a  figure  for  the  error  in  the  NCO  has  been  found,  this  assumption 
requires  eorreetion.  If  the  error  in  the  NCO  signal  is  now  added  and  the  whole  error  analysis 
repeated  from  the  PD,  the  error  output  after  the  seeond  iteration  is, 

C,ror,2  =4.5743£-4  (4.15). 

This  error  analysis  shows  that  as  the  feedbaek  system  eontinues  to  operate,  the 
error  ean  grow.  At  this  point  an  idea  of  how  fast  errors  ean  propagate  and  grow  in  the  feedbaek 
system  is  determined.  By  defining. 


C,i  C,i-i  ’ 
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(4.16) 


the  for  z  =  1  is  =  9.2296E  -5  .  This  does  not  mean  that  the  error  will  grow  at  this  rate, 
because  errors  can  cancel  themselves  out.  However  this  is  the  maximum  rate  at  which  they  could 
be  expected  to  grow. 

Figure  15  shows  the  output  of  the  fixed-point  NCO  subtracted  from  the  output  of 
the  corresponding  floating-point  NCO.  This  difference  represents  the  error  in  the  fixed-point 
implementation.  The  error  can  be  seen  to  oscillate,  with  a  maximum  value  of  around  2AE-3  . 
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Figure  15.  Error  Plot  of  the  Output  of  the  Fixed-point  NCO  versus  Time 


How  does  this  affect  the  performance  of  the  PEL?  Eor  the  SNR;  =  15dB  used  for 

this  simulation,  it  essentially  does  not.  Using  a  maximum  error  of  2.4E'  -  3 ,  and  the  input  SNR  of 
15  dB,  the  error  reduces  this  by  0.164  dB.  The  new  .  This  is  not  a  problem  for  the  PEE 
considered.  The  PEE  locks  consistently  and  as  shown  in  the  sections  above,  within  expected 
theoretical  limits.  The  quantization  error  of  a  fixed-point  system  could  be  a  concern  if 
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implemented  in  a  low  SNR  system,  but  for  communication  purposes,  with  SNRs  over  10  dB,  the 
reduction  will  not  impact  the  performance  of  the  PLL. 


b.  The  8-bit  PLL 

An  8-bit  PLL  was  attempted,  but  as  can  be  surmised,  the  8  bits  did  not  give  a  large 
enough  range  or  a  high  enough  resolution.  When  one  parameter  was  fit  to  the  dynamics  of  the 
PLL,  the  other  was  too  small,  causing  critical  errors.  This  last  section  examines  why  this  was  the 
case. 

The  only  component  of  the  PLL  that  could  not  be  made  to  fit  an  8-bit 
representation  is  the  loop  filter.  This  is  because  in  this  component  the  largest  ranges  are 
encountered.  Recall  from  Equation  (4.5),  the  output  of  the  loop  fdter  is  described  by 

Uj (n)  =  Uj^(n-\)  +  (n) -I- b^u^ (n-\)  , 

where  b^  =  635.1475  and  hj  =  -621.4397 . 

The  initial  cycle  of  the  PLL  sets  m^(0)  =  m^(0)  =  0 .  This  means 
Uj{\)  =  635.1475m^(1)  .  Assuming  the  worst  possible  phase  offset  between  the  reference  and  the 
signal  sent  from  the  NCO  of  ^/2,  u^{l)  =  7t / 2 -1.5 .  From  Equation  (4.5), 

Uj  (1)  =  635. 1475  *1.5-  982 .  If  a  phase  offset  of  -TT /  2  is  used,  (1)  -  -982 ,  thus  at  a  minimum 

this  range  must  be  included.  To  achieve  this  large  range  with  8-bits  requires  a  scaling  of  S  =  2^ . 
The  corresponding  precision  for  this  scaling  is  8.  The  second  iteration  comes  along  and,  because 
a  correction  has  been  made  to  the  NCO  output,  the  phase  error  is  slightly  less,  around  1.4. 
Plugging  in  the  formula  for  Uj  , 

(2)  =  982 +  933 -932  =  939. 

ft  is  important  to  notice  that  the  last  two  terms  of  this  sum  add  to  1,  but  with  a 
precision  of  8,  the  result  is  zero,  thus  m^(2)  =  982  .  No  correction  from  the  second  iteration  is 

made.  The  PEE  cannot  lock  as  the  corrections  to  the  phase  of  the  NCO  frequency  cannot  be 
realized  due  to  the  large  range,  coupled  with  the  need  for  high  precision. 
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Numerous  methods  were  tested  to  correet  this  error,  but  with  little  sueeess.  The 
sum  was  broken  down  into  numerous  sums,  so  that  a  higher  precision  could  be  obtained,  but  the 
initial  loop  made  this  approach  unrealizable.  A  second  more  promising  approach  was  to  use  two 
8-bit  numbers  in  the  loop  filter,  similar  to  a  double  floating-point  number.  Unfortunately, 
Simulink  does  not  make  it  easy  to  implement  this  type  of  representation  and  the  resulting  circuitry 
and  eventual  hardware  implementation  complicated  the  system  more  than  using  16-bit  circuitry. 

c.  Effects  of  Changing  PEL  Parameters 

If  not  already  realized,  it  should  be  mentioned  that  the  scaling  of  the  fixed-point 
numbers  is  a  crucial  step  in  developing  a  fixed-point  model.  The  range  of  the  numbers  and 
precision  must  be  considered  and  evaluated  at  each  step,  and  when  the  two  conflict,  priorities 
must  be  set  to  determine  a  solution. 

It  must  be  realized  that  changing  parameters  can  result  in  changing  the  range  or 
precision  of  certain  points  in  the  model.  This  is  especially  true  when  changing  the  sampling 
frequency  of  the  system.  The  sampling  frequency  can  have  a  large  range  in  itself  and  thus  a 
dramatic  effect  on  the  PLL  operation.  For  the  model  considered  here,  a  sampling  frequency  of  20 
kHz  was  used.  Reducing  the  sampling  frequency  to  14  kHz,  the  model  reluctantly  locks,  but 
jumps  are  seen  at  the  output,  indicating  unlock  at  certain  places  due  to  out  of  range  errors.  Drop 
the  sampling  frequency  down  to  12  kHz  and  the  model  is  nearly  useless  as  a  PLL.  The  PLL  does 
seem  to  lock,  but  at  the  steady  state,  is  unlocked  as  often  as  locked.  The  output  signal  of  the 
system  is  extremely  jumpy,  and  useless  in  a  communications  application. 

Changing  the  lock  range  of  the  system  can  also  have  dramatic  effects  on  the 
system.  Once  the  lock  range  is  given,  the  natural  frequency  is  calculated  from  the  given  lock 
range.  This  value  is  used  to  calculate  nearly  all  of  the  parameters  in  the  PLL.  Changing  the  lock 
frequency  can  cause  an  avalanche  of  changing  values,  resulting  in  a  different  scaling  requirement 
for  the  fixed-point  representation  of  values. 

E,  CONVERSION  TO  HARDWARE 

To  complete  the  analysis  of  the  fixed-point  PLL  simulation,  converting  this  model  to  a 
hardware  implementation,  either  on  a  DSP  or  FPGA  chip  needs  to  be  mentioned.  Mathworks  has 
a  product  that  works  in  conjunction  with  Simulink  to  produce  C/C++  code  from  Simulink  models. 
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It  provides  a  debugging  feature  and  eode  optimization  paekage.  This  software  paekage  is  ealled 
Real-Time  Workshop.  At  the  time  of  this  writing,  the  software  paekage  was  not  available  to  the 
author,  thus  testing  the  eonversion  to  eode  was  not  possible.  The  literature  suggests  that  this 
paekage  will  translate  a  Simulink  model,  whether  diserete,  eontinuous,  or  fixed-point,  into  a 
portable,  working  program. 

Onee  this  program,  in  C/C++  or  Ada,  has  been  developed,  eonverting  it  to  a  VHDL 
language  is  simply  a  matter  of  syntax  (HDL  stands  for  Hardware  Deseription  Language).  This 
type  of  language  was  adapted  as  the  need  for  a  language  that  could  be  used  by  a  program- 
controlled  machine  for  generation  of  final  hardware  was  recognized.  Several  commercial  HDL 
packages  are  available  today. 


In  this  chapter,  the  procedure  for  designing  the  fixed-point  PLL  model  was  described. 
This  was  done  by  building  a  floating-point  model  and  converting  it  to  a  fixed-point  equivalent. 
Analysis  was  done  by  comparing  the  fixed-point  simulation  results  to  those  obtained  from  the 
floating-point.  Finally  an  error  analysis  was  done  to  determine  how  large  errors  could  be 
expected  to  grow  due  to  the  fixed-point  arithmetic  quantization  error.  The  final  chapter  deals  with 
capabilities,  limitations,  and  future  research  for  a  FPGA  based  PLL. 
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V.  CONCLUSIONS 


The  objectives  of  this  research  were  to  build  a  fixed-point  PLL  and  analyze  the 
performance  of  the  model.  The  was  accomplished  by  building  a  floating-point  PLL  model  and 
converting  it  to  fixed-point  arithmetic  and  comparing  the  results  of  the  two  simulations.  This 
Chapter  discusses  the  capabilities  and  limitations  of  the  fixed-point  model,  and  concludes  with 
recommendations  for  further  research  in  this  area. 

A.  CAPABILITIES  AND  LIMITATIONS 

The  fixed-point  model  developed  was  able  to  lock  onto  a  given  reference  input  signal  with 
comparable  performance  to  the  floating-point  or  analog  equivalent.  Lock-ranges  and  lock-times 
in  the  model  were  equivalent.  The  model  proved  to  be  able  to  track  a  drifting  frequency,  lock 
onto  a  frequency  step  larger  than  the  theoretical  value,  and  a  phase  step  with  accuracy  and  within 
theoretical  limits.  As  a  fixed-point  implementation,  it  is  a  viable  option  for  use  on  a  FPGA  or 
DSP  chip,  fulfilling  the  requirement  for  a  PLL  in  a  software  radio  package.  The  reason  this  is  the 
case  is  because  the  fixed-point  operations  are  much  faster  than  the  IEEE  floating  point  standard, 
allowing  real-time  radio  in  a  downloadable  software-like  format.  The  second  major  advantage  of 
the  fixed-point  implementation  is  the  memory  storage  required.  Because  the  range  of  the 
numbers  is  known,  16-bits  can  be  used  and  scaled  appropriately  to  represent  each  number  as 
opposed  to  32-bits  or  64-bits  of  floating  point.  This  reduces  the  amount  of  memory  required  by 
50%  or  66%,  a  valuable  commodity  on  a  small  DSP  chip. 

The  cost  of  fast  operation  and  memory  savings  is  a  trade-off  in  precision  and  range.  Using 
only  16-bits  means  the  range  of  numbers  or  the  precision  much  be  sacrificed.  The  limitation  of 
the  fixed-point  model  is  the  parameter  specific  operation  of  a  given  model.  The  point  behind 
using  16-bits  and  a  fixed-point  operation  is  that  the  range  of  numbers  is  known,  and  the  window 
of  numbers  that  16-bits  can  represent  can  slide  to  accommodate  that  range.  The  other  side  of  this 
is  that,  if  parameters  are  changed  and  that  window  requiring  representation  changes,  the  fixed- 
point  model  could  fail.  To  have  a  completely  portable  model,  the  scaling  of  the  fixed-point 
numbers  needs  to  be  variable.  By  allowing  the  scaling  to  be  variable,  any  change  in  parameters 


57 


can  be  accounted  for  by  changing  the  scaling,  hence  sliding  the  window  of  fixed-point 
representation  to  match  the  changing  range. 

The  16-bit  fixed-point  model  introduced  error  into  the  PLL  that  otherwise  would  not  be 
there.  This  is  seen  in  the  quantization  error  in  representing  a  number.  To  speed  up  arithmetic 
operation,  fixed-point  numbers  are  represented  in  a  way  that  makes  binary  arithmetic  -  the 
arithmetic  a  processor  uses  -  much  faster.  This  representation  is  not  optimal  for  representing 
real-world  values.  As  was  shown  in  Chapter  4,  this  error  resulted  in  a  reduction  of  the  SNR  of 
around  0.16  dB.  For  radio  communication  purposes,  this  is  not  an  obstacle,  but  for  applications 
using  a  small  SNR,  limitation  on  fixed-point  representation  needs  to  be  considered. 

As  mentioned  numerous  times,  fixed-point  implementation  has  a  smaller  range  of 
representable  numbers,  given  a  certain  precision.  For  the  loop  filter  used,  a  certain  precision  is 
required  to  make  sure  small  additions  have  an  effect  on  the  output.  This  results  in  the  precision 
being  the  priority  and,  hence,  the  range  suffers.  Two  limitations  are  the  cause  of  this  limited 
range. 

First,  an  8-bit  PLL  was  not  realizable.  By  using  8-bits,  the  range  was  decreased  to  the 
point  that  numbers  were  over  50%  outside  the  representable  range.  This  produced  critical  failure 
errors  in  the  feedback  loop,  preventing  the  PLL  from  locking.  Several  solutions  were  attempted 
and  could  hold  promise,  but  a  16-bit  fixed-point  implementation  was  preferred. 

The  second  limitation  of  the  reduced  range  is  that  the  pull-in  range  of  the  system  was 
reduced.  While  a  second-order  analog  or  floating-point  PLL  with  an  active  P-I  loop  filter  should 
have  unlimited  pull-in  range,  the  fixed-point  counterpart  did  not.  This  was  due  to  the  range 
restrictions  imposed  by  the  fixed-point  representation.  Because  only  a  limited  range  could  be 
represented,  when  a  frequency  outside  that  range  was  encountered,  the  value  could  not  be 
approximated  in  the  fixed-point  scaling.  This  inevitably  prevented  the  system  from  locking. 


B,  RECOMMENDATIONS  FOR  FURTHER  RESEARCH 

The  main  area  for  further  research  is  to  actually  build  a  phase  lock  loop  on  an  FPGA  or 
DSP  chip.  This  step  will  have  two  steps.  First,  the  PLL  model  must  be  converted  to  useable 
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code.  The  Mathworks  Real  Time  Workshop  is  ideal  for  this  purpose.  Researeh  into  this  program 
and  the  how  eode  is  produeed  from  a  Simulink  model  needs  to  be  done.  This  researeh  would 
require  understanding  how  different  parameters  in  Real  Time  Workshop  affeet  the  generated 
eode,  and  how  to  optimize  the  final  program. 

The  seeond  step  would  involve  translation  to  an  HDL  .  While  there  are  several 
eommereially  available  HDL’s,  ineluding  Interaetive  Design  Language  (IDL)  from  IBM, 
Instruetion  Set  Proeessor  Speeiiieation  (ISPS),  Test  Generation  and  Simulation  (TEGAS),  Texas 
Instruments  HDL  (TI-HDL),  and  ZEUS,  ereated  at  General  Eleetrie  Corporation,  the  most 
promising  language  was  developed  in  1983.  Contraeted  by  the  DoD,  VHDL  2.0  was  released  in 
late  1983.  In  1987,  VHDL  beeame  the  IEEE  standard  for  HDL.  It  is  sufficiently  rieh  for 
designing  digital  systems  and  supports  a  hierarehieal  description  of  hardware.  This  means  it  ean 
be  used  to  deseribe  programs  from  systems  to  gates,  or  even  the  switeh  level  [7].  To  eontinue  this 
researeh,  an  understanding  of  an  HDL  and  a  conversion  from  the  present  Simulink  model  needs  to 
oeeur.  This  proeess  is  not  as  easy  as  a  simple  eonversion  from  C/C++.  Erom  the  C/C++  eode,  a 
Preliminary  Design  Deseription  (PDD)  is  typieally  generated  prior  to  VHDL  eoding.  Onee  an 
aeeeptable  PDD  is  eomplete,  VHDL  eoding  ean  eommenee  [3]. 

Analysis  on  what  type  of  LPGA  or  DSP  is  needed  to  satisfy  the  requirements  of  a  software 
radio  needs  to  be  done.  This  would  inelude  figures  of  merit  sueh  as  speed  and  memory.  If  the 
radio  requires  GHz  operation,  but  an  EPGA  ean  only  handle  MHz,  a  real-time  system  is  not 
feasible.  This  is  not  the  ease,  as  “IE  Sampling”  ean  bring  the  operational  frequeney  down,  but 
researeh  on  the  trade-off  between  high  speed  and  affordability  is  required.  Researeh  into  how 
mueh  memory  would  be  required  to  not  only  put  the  PEL  on  the  ehip,  but  the  other  eomponent 
parts  of  a  radio  sueh  as  a  demodulator,  deeoder,  or  delay  loek  loop  (DLL). 

C.  FINAL  COMMENTS 

The  final  goal  of  this  and  related  researeh  is  to  build  a  software  radio  on  a  DSP  or 
preferable  an  EPGA  ehip.  The  idealized  software  radio  would  allow  the  user  the  flexibility  to 
download  software  paekages  related  to  their  mission  at  that  time.  This  platform  would  be  a 
general-purpose  radio.  Instead  of  using  several  radios  for  different  situations,  modulation 
formats,  or  radio  protoeol,  a  single  unit,  eapable  of  eommunieating  with  all  formats  by  a 
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downloadable  software  package  would  be  used.  By  doing  all  the  demodulation,  and  signal 
processing  in  software,  this  goal  can  be  achieved.  A  typical  radio  signal  would  be  digitized  at  the 
antenna,  or  more  likely,  an  analog  frequency  translator  would  be  used  to  shift  the  RF  signal  to  an 
intermediate  frequency  to  reduce  the  sample  rate  requirement  of  the  A/D  converter.  With  the 
signal  digitized,  all  radio  processing  function  are  done  in  software.  An  integral  part  of  this  radio 
function  is  the  PLL.  To  realize  these  functions  on  an  FPGA  chip,  fixed-point  arithmetic  is 
required.  For  reasons  already  specified  such  as  memory  space  and  speed,  DSP  and  FPGA  chips 
use  fixed-point  representation  rather  than  floating  point.  Thus  a  fixed-point  PLL  model  was 
developed. 

The  fixed-point  PLL  model  simulated  in  this  research  performed  comparable  to  an  analog 
or  floating-point  equivalent.  The  fact  that  it  was  implemented  using  fixed-point  arithmetic  makes 
it  ideal  for  conversion  to  an  FPGA  or  DSP  chip. 

The  limitations  introduced  by  the  fixed-point  arithmetic  must  be  considered  in  designing 
the  final  product.  Understanding  the  nature  of  fixed-point  scaling  and  how  parameters  can  affect 
the  ranges  of  an  application  are  essential.  For  the  model  in  this  research,  a  16-bit  fixed-point 
representation  was  assumed.  An  8-bit  model  would  have  been  ideal,  but  the  limitation  on  range 
and  precision  of  only  8-bits  caused  a  deferment  to  16-bits.  Fortunately,  DSP  and  FPGA  chips 
come  in  16-bit  and  even  as  high  as  32-bit  on  the  commercial  market,  making  the  16-bit 
implementation  a  practical  solution. 
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APPENDIX  A 


This  Appendix  gives  the  code  for  the  MATLAB  floating-point  PLL  and  the  fixed-point 
Simulink  model  used  for  simulation. 

Floating  Point  MATLAB  code 

For  ease  of  reference,  the  variables  used  in  the  code  will  be  matched  to  the  notation  used 
in  the  text. 

•  sigi  =  in-phase  component  of  w, . 

•  sigQ  =  quadrature-phase  component  of  Wj . 

•  zeta  =  ^  . 

•  Wl=  co^,  Wn=  co^. 

•  taul  =  fj,  tau2  =  ^2 . 

•  phi  =  02  • 

•  theta  =  9^ . 

•  FI  =  output  frequency  of  U2 . 

The  following  is  the  code  generated  for  PLL  model. 
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9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9- 

ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 


o, 

o 

g, 

o 

o, 

o 

o, 

o 

o, 

o 

o, 

o 

g, 

o 

g, 

o 


ADPLL2  (All-Digital  Phase  Lock  Loop) 

Implements  a  phase  lock  loop  assuming  a  complex  in-phase  and 
quadrature-phase  input  signal.  Adapted  from  Reference  [1] . 

Michael  Johannes 
12  September,  2002 


g, 

o 


g, 

o 


9-9-9-S-9-S-9-9-9-9-9-S-9-9-9-S-9-9-2-S-9-9-9-S-9-9-2-S-2-S-2-S-9-9-2-S-9-9-2-9-9-9-2-9-2-S-9-S-2-S-9-S-9-9-9-9-9-9-2-S-9-9-9-S-9-S-2-S-2-9-2- 

ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 


%  Initialize  signal 

parameters 

Fc  = 

2100; 

g, 

o 

Input  Reference  signal  frequency 

UlO  = 

1; 

g, 

o 

Input  Reference  signal  amplitude 

U2  0  = 

1; 

g, 

o 

Output  NCO  signal  amplitude 

FO  = 

2000; 

g, 

o 

Center  Frequency  of  NCO 

Fs  = 

40000 

g, 

o 

Sampling  freqency 

Ts  = 

1/Fs; 

g, 

o 

Sampling  duration 

%  Construction  of 

Reference  signal  sigi  and  sigQ 

SNRdb 

=  15; 

g, 

o 

Input  Signal  to  Noise  ratio  in  dB 

BW  = 

1000; 

g, 

o 

Prefilter  bandwidth 

pod  = 

rand*3  60 ; 

g, 

o 

Randomize  phase  offset  in  degrees 

por  = 

pod*pi/180; 

g, 

o 

Convert  random  phase  offset  to  rad/sec 

t  =  0:Ts:50/Fc;  %  Vector  of  time  samples 

%  Calculate  noise  power  from  SNR  and  prefilter  bandwidth 

sigma  =  sqrt ( . 5/ (10^ (SNRdb/10) ) /BW) 
sigma  =  sqrt ( . 5/10^ (SNRdb/10 ) /BW) ; 

%  Signal  construction  with  noise  added 

sigi  =  U10*cos (2*pi*Fc*t+pi) ; 
sigQ  =  U10*sin (2*pi*Fc*t+pi) ; 
sigI  =  sigl+sigma*randn ( size (t ) ) ; 
sigQ  =  sigQ+sigma*randn ( size (t ) ) ; 

%  Calculation  of  PLL  parameters 


=  -6; 

"6 

Damping  factor 

200*pi; 

g, 

o 

Lock  range  defined  according 

to  expected 

g, 

o 

input  characteristics 

Wl/2/zeta 

g, 

o 

Natural  frequency  calculated 

from  lock  range 

1; 

g, 

o 

Phase  detector  gain 

1; 

g, 

o 

Loop  filter  gain 

%  Loop  filter  constants  and  lead-lag  gain  constants  calculated 
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taul 

=  Ko 

*Kd/ (Wn^2 

tau2 

=  2* 

zeta/Wn; 

al  = 

-1; 

bO  = 

(Ts/ 

(2*taul) ) 

bl  = 

(Ts/ 

(2*taul) ) 

%  In 

itialization 

phi  = 

0; 

theta 

=  0 

r 

Uf=0; 

theta 

I  = 

0; 

thetaQ  = 

0; 

FI  = 

2*pi 

*F0 

%  General 

case  pha 

for  n 

=  2 

: length ( s 

%  Lead  constant  gain 
%  Lag  constant  gain 


Output  phase  of  the  NCO 

Phase  error  of  the  two  signals 

Output  signal  of  the  loop  filter 

Error  of  in-phase  signal 

Error  in  quadrature-phase  signal 

Output  frequency  of  NCO 


phi  (n)  =  phi (n-1 ) +F1  (n-1 ) *Ts;  %  Phase  error  from  current  cycle 

%  added  to  the  phase  of  the 

%  output  signal 

%  Subtract  2*pi  whenever  the  phase  is  larger  than  pi  to  keep  the 
%  phase  bounded 

if  phi  (n)  >  pi 

phi (n)  =  phi  (n)  -  2*pi; 

end 


%  Implementation  of  phase  detector 

thetal  (n)  =  sigi (n) *cos (phi (n)  ) +sigQ (n) *sin (phi (n)  )  ; 
thetaQ  (n)  =  sigI (n) *sin (phi (n)  ) -sigQ (n) *cos (phi (n)  )  ; 
theta (n)  =  atan (-thetaQ (n) /thetal  (n)  ) ; 

%  Loop  filter  equation. 

%  Implements  the  transfer  function  l+s*taul 

o, 

o 

%  s*tau2 


Uf (n)  =  -al*Uf (n-1 ) +bO*theta (n) +bl*theta  (n-1 ) ; 

%  Implements  the  equation  for  the  NCO  Wo+Uf (n) 
Fl(n)  =  2*pi*F0+Uf  (n) ; 

end  %  for  loop 
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Fixed-point  Simulink  Model 

The  Simulink  model  was  broken  down  into  its  eomponent  parts:  Phase  Detector,  Loop 
Filter,  and  Numerically  Controlled  Oscillator.  Figure  A1  shows  all  components  in  the  PLL 
system.  Figures  A2-A4  show  the  Simulink  block  structure  for  each  component.  The  scaling  of 
the  fixed-point  representation  for  each  block  is  shown  next  to  that  block  at  the  lowest  level. 


Figure  A1 .  Phase-Lock  Loop  Simulink  Model 

The  PD  takes  three  inputs.  The  inputs  Ini  and  In2  are  the  in-phase  and  quadrature-phase 
reference  signals,  depicted  by  the  [t',  sigl']  and  [t',  sigQ'] ,  respectively.  These  signals  are 

generated  in  the  MATLAB  workspace,  and  imported  into  Simulink  using  the  From  Workspace 
blocks.  The  t  input  of  that  block  simply  tells  the  model  it  is  a  sampled  signal  and  the  duration 
between  samples.  The  input  In2  is  the  phase  error  feedback  from  the  NCO. 

The  FixPtGateway  Out  blocks  take  the  fixed-point  value  and  convert  it  back  to  a  double 
floating-point  representation  so  that  it  can  be  graphed  versus  time  in  the  Scope  block. 

Figure  A2  is  the  phase  detector  block  structure  embedded  in  Figure  Al.  The  three  input 
lines  are  seen  at  the  left  of  the  figure.  Inputs  Ini  and  In3  are  the  sampled  signal  and  immediately 
go  into  a  FixPtGateway  In  block  that  converts  the  floating-point  values  to  a  fixed-point 
representation.  The  final  input,  In2,  is  the  total  phase  of  the  NCO.  It  serves  as  input  to  a  block 
labeled  sin/cos.  This  block  takes  the  phase  and  by  use  of  a  look-up  table,  outputs  the  sine  and 
cosine  of  the  input  phase.  The  output,  Outl,  of  the  PD  is  the  phase  error  of  the  system  plus  higher 
frequency  terms. 
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Figure  A2.  Phase  Deteetor  Simulink  Model 


The  next  bloek  from  Figure  A1  is  the  loop  filter.  This  bloek  is  shown  in  detail  in  Figure 
A3.  The  constants  635.146  and  -621.44  are  the  lead-lag  constants  andbj.  The  1/z  blocks 

signify  a  unit  time  delay.  This  block  saves  the  output  from  the  last  sample  time  and  outputs  it  at 
the  current  time  sample.  The  output  of  this  component  is  the  phase  error  of  the  system. 
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Figure  A3.  Simulink  Active  P-I  Loop  Filter 

The  final  component  block  of  Figure  A1  is  the  NCO.  This  subsystem  is  shown  in  Figure 
A4.  It  takes  as  input  the  error  from  the  loop  filter  and  calculates  the  total  phase  error  of  the 
system.  This  phase  error  is  fed  back  to  the  phase  detector,  and  the  process  repeats  until  the  error 
is  forced  to  zero.  The  complicated  circuitry  at  the  end  of  the  NCO  is  a  comparator  and  switch 
combination  that  compares  the  phase  output  io  n .  If  the  output  is  less  than  n ,  the  switch  outputs 
its  top  input,  which  is  just  the  signal.  If  the  phase  is  greater  than  Tt ,  the  switch  outputs  the  bottom 
input,  which  is  the  signal  minus  iTt .  This  circuitry  keeps  the  phase  from  growing  unbounded  or 
out  of  range  of  the  fixed-point  model. 
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Figure  A4.  Simulink  Numerically  Controlled  Oscillator 
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